


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991 


Approximate confidence limit procedures for 
complex systems. 


Yee, Kah-Chee 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/27995 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 































' 
' 
. i ' 1 
1 ' 
‘ is 
' . oon 
' ', 
' i ' 1 
Fi ; 
ie hg ceed 
' 4 . Dec Fa | 
; mer 
1 ' 
eee " 
. ' ' 
: : 
1‘ om ‘ LJ 


a ee | . 
‘ . oe + 
. . 
. . “fat e¢ 
v ‘ ot = 
= * ou a6 
se . ‘ 
to. 1. 
. oe ‘ 
. 
ower aa, 
. *, . rie 
Ce ee 
= cer ' 
Mo 6 . 
a 1 = 1 
oe . ' 
« . Se 
A 1 
2 . 
on «oy on 0% 
" 1 . 
' , . ' 
‘ $ . 
' U ee . 
‘ oe ae 
Lan eat 
s : ' oe 
‘ 1 ‘ ' 
. . oo. ' ‘ 
pee ‘“ r 
se te 
se . . 
1 . ‘ . 
bs te 
ou '. ' 
. . '% ' 
. . » 
. . 
' “8 
' 1 . . 
so) ee ie 
1 . ' . 
ray . » . 
eo. » 7 1 ° 
. . q 
eo Wu ' 1 
' 
Pe . . + 
° . 
on 
7 8 ' 
‘ oa ne 
“ ' 
1 o eo «= 
- . 
. . >. = 
Ly .* Coo, . 
ve ‘ - 
ree Ta) a 
eee > of 
id é 
. ve 
° - 
ee x ea SE 
: é 
oe op 
» 
eee 1‘ 
ee 
. » 
a é ie ry 
of te . 
' t 
. Le 
. . 
* » 6 a | 
. « 
. > 44 
on "6 ° 
oe a) r 
] 
» - oe 
» . ot 
ot 6 spre "4 
ss . . ® e 7 
« toe ° 
e 
* . . . 
‘ . s 
. os ‘ 
. - 
. @e - 
> . 
Ps . 
». e 
oe 
¢ » “ La eee 
> 
bel s. 
. ° * ud 
» . » 
. é . ' 
. . 
*4 
. | ad 
. . ‘ 
F . 
' 
. 
» A ‘ 
« 
a. . 
. ‘ 
' ) ay 
¢ 
. ' 
. rf ee . 
. - 
‘ . . 


ute on ne 
Bt Nas nigh ties a of A Rite WAS rei ps sss ebay ned 
: a] | 
WD MO Be hae Ob eet ar igbateetile ar tb 
ma "we Meirson. ate 
von ia Pate 
é, 
ey me 1a 
. se wee, Elites ek i 
" ‘ : 4 at sys tee nd tiesee 9 sPq stipe 
er} i : oh OF paaeuay A Si AP 1,4 A618, AO te ae ney! 4 aot a 
. Wer ® oo4 At ow wh ae ’ ttt F ae tsotaee r ale OAs i iRaten 3% 1 3 
1 ' en rar a ‘ ae ee re a ' "ie ° Shee PEN A ital Finger Realty Sitbeaye 
CRS nn oy en “ in On @ abs at Eh I lene 1 Otte RUE Os und t uh a Vode city Bean ya 
an . ‘ ' oot og hie i : tee tsor we tet Mem ify a ¢ tpgn ak, qa! x #8, z ara U ads 
ry i ea oo "a fey 1 ne pee ify 1,9 (bitte ' eth one eget af Das ow, ee gems pu reat 
i Cora Ppt ' 1 owt the pte ‘ aq Yao + | og Re 
' hie 3 a ee one Cy, ast rsh oye ot top t ia iy 
to 1 ' 4 meg Peer Oe eee - st .f 04 25 he as 1) ’ Aa it 
F, A ' tte 2% 2 5 A ' i: . ' aie Fino a oat ‘4 ‘i mepei nas ' agl ” oN ar iia 
. A Kay an ae 8 eur toy ® a, § 0 2 08Qe tee OD, ‘ a i) it tides ire iatntate 
1 ' tate ot ont oe Lise 4 « 1 tbrcqn * er hae 1 Ua ad ' w' 1 ita a7 
: , oe ay Oey ed: ‘ ' i és U Ade? re.4 aybe te 
Se ese i aa eo eee SNe wk ihaiaes ari Jae ae Sein 
. nt 1 ay . 
. on . ue ie - va, mn ee 4 Ty Nay ' ek He aaa 
a ~ 4 + sha etih ge ham a ns age i 
' ret poet ' . 1 ow whet alge owe a qs 1°U oe D raga un iat , 
1 "arp ed prety tl 1 ' . dadedet, a0 3D mrt o ty toe "S, eid mate Ris. rye 
oye ae es ' aunt arse at ise btm oo HE ets AVN oth 
' ' HaSin ty a e SoA Savile ab 
Se ome H Ae antnn cl — 
. ¢ a 1 ale Podite cas esis ated ‘4 
a eeecity ore oe corny be bib Me CN Re dine Deny ts rea “ha Agtiees sass ne” Ran estat > 
a ae - * ats a Abed, mT law, Paap hgei a it hai cat aimee eee 
soo4 ue, ft aue ‘ he adh nea. ' Peasheolue Line jonny aris Sehy Hie sale 
' ‘ f ra ; 1 1 site Be 8 afaset ‘ veal rea iad ii tA ora mn M 
Cre be eg me Ne Sia a ames ans LOE Pg mel ead dads isa ae On 
1 1 ot 1 poe on et ec tgatyy! oenoy PP ais bee! teat - 
i ' 4 Pear eo 4 os! ae a: wy oon ge ow, 334 He a Re ten re! Z a4 sa 
Pt Aor + on . 11 bP ay tet Taree 8 99 age kh 
a ‘ ont ue 11 0 Mert 1 1g 8 Ros 
Bl siege eit % 7 vie Dit gt erg on alk Belen fart A. rife 
. 2 ' ' 1 1 ed oa ' e? 190” ree! 
rhe or ae oo te ye U Ue Ley ow ' .. age hee ' +o Senge Pie he : Vag ee 
ite eT i sas a) son PF Wy he el eine * $8, arti 
' alae at dred o fey 4s 
.imrieake nee 
te rar Tee se Se onpaesie a tholecsa iY 
: oa, be ‘ ct I aby P . mts be - 
ia eee Hi hy 1 og'’ Tia mit: Bek oi brh Scejee % Sees auc 






















































i rs © .Gittsd 4, on Teed Sani 
a oF ae tee toe) car er) i We | toy ae : Vues ors at Desay: vasanys + ow ery 
“ae wig, A) es By Uae TG ee eb wangett * . ce ee At ba 3 yeh i ie Hy 4S eh pre 
e J 1 wae . : 
. ght eB ite a. sie Oh te sf en a ies Fp nat ug 
fa f Saf 
oe de ' Aut Ae Cer : i ped of Bie ‘ aa z 
nat MU ; Sai 2 , Ue ast ‘mes, ‘ me ea eC} Ly sarees 
as ENS : eae ae eee Ngee pad eho ob Me et end ewan Una d P an en 
‘ ' 9 - t weoat 4° ' a 4 t 





a ea as Ory amet “° ge 
» 4 ta a : hae i ayer 
re 





ening at 


oes 











CW etyaret 
Part eh 
































cra) a 








a Q 
i . ree f Pee eS Trae Py Soe 3 atl 
oi - LJ «! 
rte PAY! etesegihs # As bu 







































































































































































































































































































of My 


e 
a HSN Tae ved hea 


Catt 





AGost 





Sree af 
oh ing wae 7 ia) ‘ a. iv mt a ' alte 
eh oder r : os aa soit yy) 












Shirase 
hea aet tae 
- YA hese 

r ar 
Rabe 





ite at 
' ' 
ootyee 




























“ap 2,0, UU te dorky, 
yar byandt ttl Ay, 
Whe aptag Rh 
















































Wat 4 J a Pol 
Hey! ‘ ‘ se * = r ry 
Ga ms an 2 aie i, unt snat sion Sour ree 
5 ea ote %e ae a a 4A ae saa her fab te | ; + Saree ke by ee 
me a wits fe pe Pied euhte fan's ff Nee Aver acct) 
ny , “bent 









bo mae palbats! 


Yat sae 
i eee Fostelet Hey, 


Ae 





: pine : re 
‘ . i 5 r 3S is " ae ae 
Fa AS goseg Pipe ese ce etek - ' kK gy Pe 
eee en ee Rew ‘ 


ef . oe a 


































































og 
a 

Petts " geht as sath 

ie - ise 1 a 5 % 78 Me, e Festhalle vate dS rr Ay 


: », yee Pe es 
i | gd, te 


Ate 


jo ALN oy 
Lit tn ¥ dite 
a0 pte aria 4, whe ped s'wrne 

mo *aidge Ma 


tte yt iY 















se sti 
ao ro ee 
ae hee 











ryt, tay " 
by tara eth of ean « 
tyey tt sous oe ip, Pasta) ee 

OF Aitee « i Se 
. 










a rris 5 
af hes 







Pe ee 


Se vest ay 


“a a tah. e& 


.” neta * 
















sebult att ott, 


wee 
t aay Pat Manone 



















ent sptat edi kha 4 
Tye a oo aut ron a ht ¢ 
int 58 %.. 
whet ee - Som ge 4. he is todas 
baste ch a ef 
wi ., rad 
Theda oe i ae, 


*s ore joe? s 
ah to he, ian) 











el PSL oe 


regi 2"st rege, 
fee 5 we § 


f eiF 7) rT 
nighapa #3 or el eeesy ey 
arate og cra cweyreateeesy 
. Me ar en, Ee td Cieiha 























































































































































































































































































































































































































































































































































eet ket 
. : . +t ++ Z ROG ee Baia. . a ie net 
' . 1 . . Pa e8 Pe =| . nal 
ning » i @ Ba ais . . as ear ti F * ‘ ee eet ee, rybpt eed 
. . pan t . wot . . ’ . ' < < 
Pere A tet ee Ger A Pater eh " a 1 1 chaeare ate "af ware Ta rae a 3 i sare oreo te ie 
7 - “eae :", ., ‘ IS Raa rcs eae CS EE aos = igs ¢ re 
. oa Ree @ ' 1 . . i ” pee is ie 9 
: mn ' fae ee j a, : aor . woe coy : "e Neat 4 ate nA rer wie olny is i oe cast 
‘ ' . ee . .- 3, fa 
. % $ . “-¢ : ee ti ee ee inst Nn Sr : ' ‘ 19 si oma, Mt R! v fag eh etetagorys ony ie lagdese aie wate it ftente 
* ou eas oer ¢@ a Hy ten Pegs ogee Be la aa 7 pe ote a sayt fie 
' r' PT fs eal eee is cuter Fra ay ae ote ea? yg" AL aaa t pre tie i ptr rena a 
‘ rT . or > " e 
oe 4 ; mae oe ee et irae 4 ri ras ae . ga a Wye y ae - ‘ if J at ‘ Reyer wea ot tbe fe “st ynpel sty" ¢ De 4 ° Aa Ay ae a tNogt BE 
» ‘ ‘ te 2 ' en a oe Ye \. ‘ ' poe 4 seyeriom ee af rae) Ea ae Se ie * a i rhe ape gt pe ae 7 gy ecg 3 
. m4 i’ oye wen Sound wo plan gine’ ir . yk PhD abe weanes gt 32 a7 $h re 
, . " tee ' sufen & ciate 1 1% 4 % cee a oe Re pe ses PLC il sigh 9 “ie : Pse'ye vgs be ‘3 vier sea Sate q a we 
‘ 1 ' . ee 1 ' . e « 1 ai Qhete t paye' Dep nerd oe ody fs webs tare ee oF it £ ey oghans = ne " . 
F . “o ie H 3 1 rt ot t& *; ' ae $y yee ri ror Uh 20 = ahibee opti otal a a peat 7 bef i ea rd ef 
. ete. Bre ', te as . i " , Sr i ' oes % y inte eH bel Pe das kets ‘ple. “ 
Rig ieke - 00 : Py y Peeeze hoe hey, haa i 
‘ RC oes ee Vaccine ge Mere a artes aan LCD nL) ee OE, dee Wigihecia te Bee a fide Wap ea ha pet ety 
. . yy v7 . ' teres ' 4 . os nn vrei. a Leda Sahed ats Pa 
te? " erm 8 Pa er) ci thee “bh of ¢ Sr wahy fshanes grey vt 
= t . “a hey eid ie .. : a0 aa ae : rf om : ccreaire im + : 274.07 rete 
. ‘ hye ¥} 
ine to Dieter! wr6 lt oa gr elt a) Sone iret Hani Yackstis eb Ard la fl bere ch Bl te oe | rtyegiges teuer Pata 
: Wresyorgh ten oty i? 
Lc pa ae . git” 37 ee aan 42 guy Bess wines ca mat er phe Jn Pepdipeyeg a roe dae seeie ge Peey e 
te oe fee ia) oe ng hoo nay ee Gewese, SF ye wipe ¢ ne etae ef yn erent ieee Eaerk ab SLL 
or . i ’ ay a) tA OR tte nee . te » 8 14 Vay reise ete roel geetery yoy epee ath eer ercr’ I ed 
Cite mer * . ‘ of mee : © « 8 gas 14 4 bo yip wer te By pera eat ty bs ees > A 
ob an : ' id Nata ha a? ' eg e lite ee dele, 4 
: a . sweat eet ‘ . . te ae ' os ? ' eae hs UA et => 1 "of ‘3 ve “"Y 
. . . ef 
. a) 8 . oh. ae t ee | ee he | z ’ re aL) a v) ‘ the wey F eRe Sac 
' vy ' wate a » nga Carne iL miige s WH one ee *e "a> “ol ve rel ey ” inp woe on a ! yeaa’ yc Lea op then, at Ay e's i: Ste Feed gL rears Pad eananue neaeelee € 
on . . 4 f Yop awe’ ow Oo wee Ore Wg Pes lige gto f P ¢ oR ae Pw ee ela e 2a eres ” ets : shetenatatyth stated estat te ery" ite se te dme poop oagr ged" 
Le bod es U . oa, See SS 1 oN Lal | A eee ehee : $s “pn? boas fap ere e a Be Sori tere fe oh 
’ . =. >. 6 "ae teotne =e Yr a PC eda | ol an ane “ey Vi we 5 ee rr = ind o taeneet weir e 
F . We . . voy a ’ , ar slog Cornqem’ Gh ide iat a 2 
La eu 7 ’) t's a 1 ‘ee . : DP puts arene ‘y bd . netiat 7 et Perera Pee if a loan 76 wr en ly: 5 id tf ° ti dale easy Deg de bel) ryt Eset ar erate re eat 
. ies ike ere edclae ¢ trou ene DD EGA rede nema” ys Pg? : 0k ‘iy aces erat fete ee es ote elgpateee tates 
wha str : oe gn pier eye a paeter ae meee ST aL Pee ee Mag yo Eesehe "ey Pabst to rattles ty! gee Sidastomet hart oye 
2 Weaeeno e A Lee te eS orc " Nees a arelyt tesa £0: “gp, ory ° ag ol, ye As ete: apa Paani taser aa ae gee fame he gr nly Sorooae wedst 
. Pian rye nt & . va gig OB ip bats Ve ae 4 LS ren see otsts ai 
m1 3 7 ‘ i. ice Pigs eansi ee ! + rd hae vi werd ey save yurue y sor Morr oes Hs . ed 4 ee td Sate 
' ’ . 
tS >t oF peng (ene 5 iy von ey 8 
8 1 te. . ' Se Ce | ‘, ¢ z vy 77 pes ee" A *y pe Lf te Taint pee hh at 
ae G keke nite! epee! f tee are BD gt PD Ln Fae Set widbinaemin Bote yensprg” oO MG rbrvegeg 9 oe 2a? 6 rer paiee- “pune eG rz. 
' . . an ‘ ei oun a tly oA, oy ter ryty TE RE es Sp.) vee At 
2 . ie ear vat AC | , 23 tos e i o'g els » Hj Ey dene eb rae yasnowe wriagy a th ek bad why! ore rhe aes oh oles iter yn ees 
: f ‘ EPATRO UN & ot, bite te =z ao APL vice m4 dpe eta Loeb! ere bees te be aety ie tir a ener hte vee 
‘oo we Sos 0 nee = (sia 8 ot ot, is i oie or" ol a Ee) ve eae $_ v<msh ots a> pe eexs aeyee’ 
' preeta # we " bs AM * OMI Ty 
’ B is ps 9 yea Fie i i ae eae a ke id Bes rere le cr 
F rein * ; car — ve oe pone ie rar y i! wie’, ' ini ire oe of m8 Fs ‘3 roy gpecehh ts ot a ey eae Are eters ree Wein Wet be, meertee ote vad Sst cay 
. Pe tp Pon te . tal Pee ye ee ee ee '‘ Ory ele o ' “1 oa H ipede sy : tetris a i Ba vere vtennte- be Pid oF a, out he Ie ta le a 
5 " siete 19 aie Le ofa tee eet Nears Lelatye 5 rave a ny ats ry 5 m gre eae pus ey, Pe tts Meee Seen ear Le eweyietete 4 4 +4 La 
i Ces Boe ' ° te a ie ‘ ° A ee ‘ % : he) a ate AY oy ue MY Eat sae cett ete ele oe ebaen 3 Aw) bed fi I (opricncert tren Leta ace ty 
. : se » ‘5 a) at rg tee te ’ » vow Saat a tyiot oP me Ee Meh ey nipe SN athe oe Mike 
. ' v Cie at tn eee eo, s re) ef . ” a 4 tage eu ® sag na teu! * ROE Wl pw? er Mpsskete a Se Asha ta 
o> ee or ry. ' ' MP ewe Meser se s ' feoye Satiee ts su we r sik oI TADS Fixpeisemsacetece Alkereh ol 
oe * i) U ey oD Uae tk Othe 4% a bl WE Sale as tae} Teh pola cucate | ity you ey Kteobsete Nee cane 
as Satay ° Diese ge Pe eel oe CS fe Seeoeemees alee: Rigas pigilasin gute sata catae Seven cs on stistele ate ste: “ 
. . ‘ . emmy Py <, rc be 
ee Deny ie fa ' fie eee Py a P ce 1 ota tie a ie 1 At, SPimegre Pi tacacangeeen en Sa ie dod a eh Bg C 
» ° » . . “ ' 9 Pye a Fi ’ 4 eee aairiy) wet i ff, keys op pte ye cB 
Sate gh cee ene er aed Gemeente acre ee i, Ae eae eee sas y ae Wriey tb $y 98 tat putonwraconm tc rcyaned weber b spre tevte! 
>> * eo teas - > 8 ip feats eV bee oh pa om $0 tty on te ¥fa ra i a w gre ‘atgenels Aa i altetytitetarty Diceaate wit ySpratamhaee sed 
Pe iT Rae ea a CL OC ‘ hsv Se 1909 de twos vie Fata 1s a NaC rere Pape wate ny ba ede he ene ving ou panne 
' te a ; ay » wait . Rs ® Fe a me a o6 aa? Acliy Lol AOE Ae a Ura Lao spel fa bk ato err mer yet te nents rae o"ire-¥ Re 
baad ‘eeren a aie ‘ vs a a : ' ie "sides 6s o peers at LR wg wap heoeees eats phe wis wy Phanus eee. cat ter awevee ss YY 
: : ee +) vf aye i e Te ek ‘ait in ’ ary ser poeeate nes weg ie Bin mig yeri aioe ronsocen tree gey ater alana ene % 
° "ee ony on r4 r aura 
Ay a ' aes cee x . e » ane Ue benan Chet in a ae eae A ° ar “hy Rene we ‘ LMG T Uda ve ww stn! Sipce Col hen ere hes be were She hk at Le 5h | dda} orn eters pipe 
F "ine oo Cea f eRe) can ones 4 Pages 5 0t Ge te a8 BS ale by rgae 089 Gary? sew Ne Oe @ aia tee LaMar Ne 84 gemeretacy es A Cabphbh tdi ie? iia tee atta fete 
of mu. te s ° ‘) . Teer 0g . gt te ay oa “teen feet ated le oF eee " 09s 0D wee Fone Lida - Beth rer ocy . o ences ardsee 4 eerie 
' . iL ame gh ies a i : Ey ba Me tall DAI a te at a rt es ii eae ~ ie! ae eden ee rg'er owe es 
‘ Phe a . bg Bi a . . : fae eC 4 2% ee ttf fas Ae y ri nye % ; a aay & eae i she 2 ah het as a pa aA Bh At 1 Hen, Ny Mpa BAA fd bet sane 
' . Aber o 4 te ‘ “by pibedo ad myetepmytonn ange eer i nagcaig & Obey, Rivets a ey Nh py aca ree ine rae pighche See ee ala LM 
* oe ' T iceent i | a . . F ot eed say th 3," aeet Fi spky eseeh i mm, 6 ee a OS A a ee ot eae le th» te eter ime e-em eh. tm <pemen 9 oe pee feng ye 
. ; ses ' oo te 94g“ Oeaty a° wae oped a'ptvee un ase! bom sot Pro ales Bye an ivtsteG tite. wat tayca rit pm" 98 gy le +g Sie ca tems soe atest o guinventsmnergeern Precise tae 
»* . a o86e {Pte tre week ‘wan a OP tence Se RS se ys whee ‘ rerary at a des oe abt bab" YPe~ Sue 20 ow" BEM ad sre 0s SA Nin Pi anere rein yas vateard Serutprures gresate'e B 
% e Oo m Nos ~ gt eee bos 4.! yi: a Ie eC ee ea Ba ee Se Pe oe (44 if ies tangs pete es 7 iasarana ye 
.? ' P oe na oo biune: e ' r) 1 a ocrchats Oa a th Ye > yee Sat APA Sy Pym pT be BVT H CT a ig Cicer ih eas tye as Lah Pier NH oe erppes. 
zai ee Dear LI ; af arn ‘ OF »1 . yh Pie eo De rai or Fe ats nt to ote AO A ae ys cay Be “we “atioath one as “4 Ra none rey pis Reade pL di wabpis ye om # s f eraeats 
« ry . . . ” 7 ' a efe's > ge “ho ' 1 f , | ™% mae 
. *! anh ees eich para Sire ttn Agr ore 1x ) done Pee a ’ Hi Le dae iter? as } 2, Wipuites Dit at had ih cwuery Sey Siete nk ta ee ae re 9 oats Sie oe eae 
» os . ° o* 8 9h © ht ON oy, C8 hee gaeziiec ee ier ote t abe oy ues alt we i a pore fv “Vesty aera *a? 7 i F Mer ‘ “3 Ad sanyo aes : See Wegn ope sa Tey 
af + de tie Aa AEH Si ye gdh de oo »% Gh 2 an a6 Bg: ae OOP Oa ee a é, ees Seta arp aati ) aye MTA ’ "ey mrad y} hectta 1% 4 arb eahe a 
. . . . Aid r) ’ *% ‘ @.* @es . ip. ilps ie wee 7° - ai bel an a. ~] “Sure fh I e 
, ' rat’ te yes, ’ = je ten te © he A yee s 1 fide *pase'} 9 ge eae jer “ele one “ecole he 
hs ta an i: *A, " : » ' ts 1° a 7 nn Ly wf DC ye ' ue tiie oP ¥ ‘saa oe ‘ ‘ te y Le, ay a oy ne 2 in Aaa Gewe yr ite : SPU Tatts ie Hacer sehde 
rie oa 4 went Meg why ~ ty eli « « + wow ey Tyepeirnta™ story re vhnat enh Ih eh fe De lt lk 4 te 
eee 4 © ape . ' 4 cS oes ¥ Foe >» Py ser pee im mbes Ft Ms Sagal gM lesa O'R ed eed ee. *@PPewr ee h pm > 
. ; » A ® n ’ ee ‘ ee as ws « IC « “he 3a B “irks FE Lp lee ah cree Rance ines sane 
« : “ee Ao . 4 as . “4 totes wand bg a 2 ! ua pa bia 
an "1 hee ' , i at . — hier I a *gaie is ia 4 %? ay, Coe cae o4 Mh teh pe Paar ae Nee rata Waist “gp &- res eine ie tip are? 
0k 2 NC ee eee pth, * Dyrs agin aie “a alate nt Rabi ena tieng- el nO sea a RRR aa a Raid an 
‘a 4 \% ™_ bed e ca . 4 
i a 5 . t 1 , 7 . % A tee ’ me Seaver ie fs phe sh od ny ‘a Ae ne “} Pater erly a" Lae a4 * re etean es Hid I f Et he sieges “9 re tik) Senet ee. O) Mi b oan ‘atic rete Po ne dipe cele 
. ‘ oe eta am nue mn 8 bey 8 ic Cheer Gana © Pe ad : aot Pe Koay We vonerer, fe nn Fd de ld lb ACA de ity ar Gat Sars “ph ie: 
ao fj - ar a J A rE z ghee ‘oare 7 a, it < gi 7 » brary yt we Ve ve Gretna fi: ao a Sfosae eh pa Past " % wean Mn a “ti aft ee? Ret oh ad ‘ont 
z . ee kOe C % 6 . ike at pais ed ihe -° NN rin uae cis Sateen 1. igh stabi Sita ath a wy i ara ho ep 59 RIS iatatye vet 
me ' : ea seers ; we 68 tae aie on wa ee ey woah apenas eo ise 1 ye Oa Tah tage betel ecpatey * cage: cy ratqed Oy enear Shen saiooe 
rad came fae ei pee = = sie meastrg . Ase eee ies Te Pettit 4 y aie 1 Ste rechent ety we & peed heat awe: Stiee sh iste Paty! oe oe ain) s cwesee 
. woe % « ‘ tee 003 tican Stare hse sip : fer eieeecea eee ton Ui aaa ad “eee ey ndiiakan 
ae cre ae ae Ps aed A oe gt >t. ae eee "Op EL ae ey oe ig Gityts hd ei te sages sea acti = Tyas yee WH ws A Pre setae gs trey sate 04. seb 
' \ 4 . ie Seals Set AD <4 oe wee ur ght Sotetee te 4 , By Sty ek a elect wi ert? Stew bod acstit te, ile veny re ahh Atay (™ she age ahs aaah ie Wels ae ane 
. ? ae AS ? . so 4% wi a or Sry ToS bets , rave dros te vie usdoe's! 7s wt} . om 8 amie micerenys yy" “tasepehy! fureeh ° 4 Seem 749 ©) es =p? wees 
Ce ee Ce aaa “ $e ae me aa Baie & "a fie ot raat gre - ot "8 sa * Fs oan ty arene eS err $7335 fer pe Saupe eet rely Khe marerceb sree eae 
: ‘ - Fe > mer 6 t a re ° a ae Piaaay a 5 aoe” 
z he bit » ¥ ry cred ° « *s f ! a5 1% ‘ os oe te ‘Eye : t aly iS aha nas fons “in Shih, anaes. NET ae 3 emyeed be rig bi 4 Seta aeth te Batol anos they 
4% os #1 . oe) so) euueas Jy he ye 7g ies ad one a ; iil i Se oe St arty y ope h rate et alate Anne ratte $ sa i typ em 5 3 ee ates’ Gumromonraey 
. fae OS t> [2 a ee ty Cc o Le, fe ie * ee ee ee aperte at ea 5 ante SM Re he ot aa) on) sae 1 tai hus leat Sg eae ba hs a 
Bp ® s t . re « . 42- « s a LE es | % 














NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


APPROXIMATE 
CONFIDENCE LIMIT PROCEDURES 
FOR COMPLEX SYSTEMS 
by 


YEE, Kah-Chee 


September, 1991 


Thesis Advisor: W. M. WOODS 





Approved for public release; distribution is unlimited. 


T?593273 





sanity f 


SECURITY CLASSIFICATION OF THIS PAGE 


REPORT DOCUMENTATION PAGE 


1a. REPORT SECURITY CLASSIFICATION ib RESTRICTIVE MARKINGS 
UNCLASSIFIED 
2a. SECURITY CLASSIFICATION AUTHORITY 3. DISTRIBUTION/AVAILABILITY OF REPORT 


1 Approved for public release; distribution is unlimited. 
2b. DECLASSIFICATION/DOWNGRADING SCHEDULE 















4. PERFORMING ORGANIZATION REPORT NUMBER(S) 5. MONITORING ORGANIZATION REPORT NUMBER(S) 









6b. OFFICE SYMBOL 
(if applicable) 
OR 


7a. NAME OF MONITORING ORGANIZATION 
Naval Postgraduate Schoo! 


6a. NAME OF PERFORMING ORGANIZATION 
Naval Postgraduate School | 








6c. ADDRESS (City, State, and ZIP Code) 
Monterey, CA 93943-5000 


7b. ADDRESS (City, State, and ZIP Code) | 
Monterey, CA 93943-5000 









8a. NAME OF FUNDING/ASPONSORING 8b. OFFICE SYMBOL 9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (If applicable) 


8c ADDRESS (City, State, and ZIP Code) 10. SOURCE OF FUNDING NUMBERS | 


Work Unit Accession 
Number 





Program Element No Project NO 












11. TITLE (include Security Classification) 
APPROXIMATE CONFIDENCE LIMIT PROCEDURES FOR COMPLEX SYSTEMS 


12. PERSONAL AUTHOR(S) YEE, Kah-Chee | 


13a. TYPE OF REPORT 13b TIME COVERED | 14. DATE OF REPORT (year, month, day) 15 PAGE COUNT 
Master’s Thesis From To 1991, September 129 


16. SUPPLEMENTARY NOTATION 


The views expressed in this thesis are those of the author and do not reflect the officia! policy or position of the Department of Defense or the U.S. 
| Government. 


17. COSATI CODES 
FIELD GROUP SUBGROUP 


19. ABSTRACT (continue on reverse if necessary and identify by block number) 








18. SUBJECT TERMS (continue on reverse if necessary and identify by block number) 


CONFIDENCE LIMIT PROCEDURES, COMPLEX SYSTEMS, LIFE TESTING, 
SERIES SYSTEMS, SYSTEM RELIABILITY RELIABILITY 







Lower confidence limit estimation procedures for the reliability of several systems are developed and their accuracies evaluated using computer 
simulation. The procedures use test data on components of the system which can have failure times with either exponential or Weibull 
distributions or both. Testing scenarios for the components can be truncated by number of failures or by planned test times. 

Although the evaluation effort was focussed on series systems in this thesis, the procedures readily apply to other systems as described in the 
thesis. The evaluations demonstrate the procedures to be quite accurate when sufficient component testing is performed. 

Two FORTRAN computer programs were written to perform the evaluation. They are annotated in Users’ Guides and can be used to determine 
the accuracy of these approximate lower confidence limit procedures for a given specific system and associated set of input parameters. 





20. DISTRIBUTION/AVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION 
EJ unccassiicovunimiteo = [J same asrerort [CL] oticusers UNCLASSIFIED 
22a. NAME OF RESPONSIBLE INDIVIDUAL 22b. TELEPHONE (Include Area code) 
W.M. Woods (408 646-2743 OR/Wo 
DD FORM 1473, 84 MAR 83 APR edition may be used until exhausted SECURITY CLASSIFICATION OF THIS PAGE 
All other editions are obsolete UNCLASSIFIED 


l 


Approved for public release; distribution is unlimited. 


Approximate Confidence Limit Procedures 


For Complex Systems 


by 


pee Kah-Chee 
Operations Analyst, Ministry of Defence, SINGAPORE 


B.Eng(Mech), National University of Singapore, 1987 


Submitted in partial fulfillment 


of the requirements for the degree of 
MASTER OF SCIENCE IN OPERATIONS RESEARCH 
from the 


NAVAL POSTGRADUATE SCHOOL 
September 1991 


ABSTRACT 


Lower confidence limit estimation procedures for the reliability of several 
systems are developed and their accuracies evaluated using computer simulation. 
The procedures use test data on components of the system which can have failure 
times with either exponential or Weibull distributions or both. Testing scenarios for 
the components can be truncated by number of failures or by planned test times. 

Although the evaluation effort was focussed on series systems in this thesis, the 
procedures readily apply to other systems as described in the thesis. The evaluations 
demonstrate the procedures to be quite accurate when sufficient component testing 
is performed. 

Two FORTRAN computer programs were written to perform the evaluation. 
They are annotated in Users’ Guides and can be used to determine the accuracy of 
these approximate lower confidence limit procedures for a given specific system and 


associated set of input parameters. 
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THESIS DISCLAIMER 


The reader is cautioned that the computer programs developed in this 
research may not have been exercised for all cases of interest. While every effort has 
been made, within the time available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered validated. Any application 


of these programs without additional verification is at the-risk of the user. 
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I. INTRODUCTION 

This thesis develops approximate lower confidence interval procedures for the 
reliability of complex systems using test data on components of the system. The 
accuracies of these procedures are also assessed using computer simulation. The 
procedures can be used for any complex system whose reliability does not decrease 
when the reliability of any one of the components is increased. 

The failure times of the continuously operating components are assumed to 
have either an exponential or Weibull distribution. Parameters in both distributions 
are assumed to be unknown. The Weibull distribution is used to model the lifetime 
probability distribution of electronic components with non-constant failure rate 
functions. It is also used to model the lifetime probability distributions of mechanical 
devices, since their failure rate functions are usually increasing with operating time. 

Lower confidence limit estimation procedures for system reliability are needed 
during the development phase of systems to provide indications of a contractor’s 
ability to meet a stated system reliability goal as development progresses and the 
results of test programs become available. These procedures are also needed to 
assess the reliability of systems that have been operating in the field for some time 
and have accumulated histories of failure data and unique configurations of modified 
or repaired components. 

Few textbooks on reliability treat the problem of system reliability interval 
estimation. Those that do usually limit the discussion to series or parallel systems. 
Moreover, the procedures they present are not adaptable to other more complex 
systems. Mann, Shafer, and Singpurwalla [Ref.1 pp 487-524] provide one of the better 
treatments of a variety of these methods in Chapter 10 of their book. This chapter 
provides an excellent summative discussion of the many procedures that were 
developed from 1954 to 1974. However, none of the procedures reviewed in their 


book can accomodate the use of test data from a mix of components with both 


exponential and Weibull failure time distributions. The procedures presented in this 
thesis does accomodate this type of system with a mixture of different component 
types. In addition, the procedures presented in this thesis can accomodate an 
additional mix of components for which only attribute data has been collected. 

Procedures developed in this thesis are extensions of a procedure developed by 
Myhre, Sanders and Rosenfeld [Ref.2}. In their paper, they assume the failure times 
of continuously operating components have exponential distributions with associated 
failure rates, \;. The test data on the remaining components, the number of observed 
failures f; inn; tests, are assumed to have Poisson distributions with associated means 
n.q;. They assume the ratios 4,/); , q,/q; and 4,/q; are known and develop confidence 
interval estimation procedures for system reliability that use this information. They 
also show that the accuracy of their procedure is not very Sensitive to moderate 
inaccuracies of these ratios. This suggests that it might be possible to estimate the 
ratios from the data as part of the interval estimation process and not suffer 
significant loss of accuracy in the interval estimates. Estimating these ratios is part 
of the procedures developed in this thesis. 

This thesis also provides an annotated computer program that can be used to 
assess the accuracies of the lower confidence limit procedures when applied to any 
specific system. Sufficient annotations are provided throughout the program in 
Appendix C. This program provides the user with a means for verifying the accuracy 
of these proposed lower confidence limit procedures for his specific system and 
testing program, that is, sample sizes and type of truncation. This capability will 


allow the user to answer many ‘what if" type of questions. 


Il. THEORY 


A. Interval Estimation Procedure for Exponential Failure Times 
A system is defined to be quasi-coherent if an increase in reliability of any one 
of its components does not cause a decrease in system reliability. The components 
of a quasi-coherent system do not need to be statistically independent. However, 
throughout this thesis, it is assumed that all components are statistically independent. 
Suppose a quasi-coherent system has k components and the distribution of the 
failure time of component i is exponential with failure rate ;. Then the system 


reliability R, can be written as a function of 4; ,7 = 1, 2, ..., k as follows: 
RAG) Gael eames 0, ,f5,.,f, ) (2.1) 


where f¢, is the operating time for component 7. Let m be any one of the k 


components and r, = d,/X,, , fori = 1, 2, ..., k. Then equation (2.1) may be viewed 


as 


Jy a Oe cr es eel os Co a (22) 


If the rs are knownand 4,, 4 were an upper 100(1-a)% confidence limit for ,., 


the corresponding lower confidence limit for R,(t) would be: 
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Specifically, if we have a series system of independent components, so that 
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If n; items of component / are tested until f; failures occur, 7, denotes the total test 
k 


time accumulated on all the n; items, and F = }°f, then the expression 
i=l 


k 
Io) Daneel 
i= 
has a Chi-square distribution with 2F degrees of freedom. See Bain and Engelhardt 
[Ref.3]. The corresponding 100(1-a2)% upper confidence limit for i,, is 
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where xX’, 7 is the 100(1-a)th percentile point of a Chi-square distribution with 2F 
degrees of freedom. 
If the testing on component i is terminated when a total test time of 7; has 


been accumulated by all n; items, then the equation for 4,, 44, becomes 


2 
§ _ X a2(1+F) 
mU(a) ~ —~-—— -ves(2aD) 
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In this case f, is random and so 1s F. 
If testing on each of the n; items of component are tested until a planned test 
time or failure, and failed items are replaced immediately, then equation (2.7) will 


be the exact expression for 4 . If failures are not replaced, then equation (2.7) 


m,U(a) 


is approximate. See Lee, Bain and Englehardt [Ref.3 pp 486-495]. Department of 


Defense document NAVSEA OD29304B "Reliability and Availability Evaluation 


Program Manual" [Ref.4 p 5-42] provides nearly exact procedures for Nes when 


testing is terminated by planned test time for each item tested and failures are not 
replaced. 

The values of the 7,’s are assumed to be unknown in this thesis. When testing 
is terminated by the number of failures, a nearly unbiased estimator for r, is 


p= Ae 28) 


where A , = G;-1)/T, and the index m denotes the component with largest value 


of 4, . The ratio, (f-1)/T;, is an unbiased estimator for \; (see Appendix A). If 1/4, 
were unbiased for 1/A,, then f, would bean unbiased estimator for r,; Replacing 4, 


with Xf /(f,,-1) in equation (2.8) will yield an unbiased estimator f, for 7;. 


t 


Multiplying by this constant f,,/(f,,-1) is nullified by a cancellation with the same 


constant in the final equation for the system reliability lower confidence limit, so 


A 


equation (2.8) is used to estimate 7; Using estimator f, for 7; , equation (2.6) 


becomes 
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It is important to note that the index m denotes the component for which 3, 


= (f-1)/T, is the largest among all the components in the system. The corresponding 
equation for the 100(1-a)% lower confidence limit on the reliability of a series 


system Is 


k 
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The corresponding lower confidence limit for the reliability of any quasi- 


coherent system is given by equation (2.3) with 7; replaced by fF, . 
In this thesis, equation (2.8) with 4, = f,/T, will also be used to estimate r, 


under exponential assumptions when testing is terminated after an accumulated test 
time is achieved (truncated), and when at least two components have at least one 
failed test item. This is done because failures will not be replaced in the time 
truncated test plans that are simulated in this thesis. In this type of testing both f, 
and 7; are random. Under time truncation, it is possible that no failures will occur 
on any component tested in which case equation (2.8) is undefined. Also, if only one 
component has one failure and the remaining components have zero failure, equation 
(2.8) would be zero for all i except the case wheni = m. 

All of the confidence limit procedures in this thesis have a common special 
method for computing the lower confidence limit of system reliability in the two 
cases of zero or one failure. This feature amounts to a modification to equations 
(2.3) and (2.10). 

When either zero failures or one failure have occured among all components, 
the test data is examined for each component to determine the total number, N, , of 
equivalent component mission tests (i = 1, 2, ..., k). These N,’s and the system 
configuration are analyzed to determine the equivalent number of mission tests, N, 
for the system that would have occurred if N,, N>, .... N, of these k components 
were assembled into systems. 

For a series system, this N will be equivalent to min { N,, N>,.., N, }. The 
100(1-a)% lower confidence limit of system reliability if zero failures occurred is then 


computed directly as follows: 


ROne = ‘le oe Ae ly 


If exactly one failure occurred among all & components, Ones will be the 


solution for p in the equation 
Pp N + Np™1(1-p) = @¢7 eee (2.12) 


These confidence limit equations are the standard binomial lower confidence limit 
equations. Equations (2.11) and (2.12) are part of the set of equations used to 


compute On. for all of the time truncated interval estimation procedures in 


this thesis. 
It is important to remember that the symbol 7; in equations (2.7), (2.8) and 


(2.9) denote total accumulated test time for component 7 ; that is 


_< Woe 
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where 7;, is the test time accumulated on the jth test item of component / and n; Is 


the number of test items of component / being tested. 


B. Interval Estimation Procedure for Weibull Failure Times 
Consider a series system with k components. Let the time to failure, X; , of 


component 7 have a Weibull distribution with density 


He) = Bt exp { -At)") , +t > 0 » (2.14) 


Then 
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and 
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where A, = .°', X., is any one of the A, ,7 =“, 2, .... kK fand7, = /A= eine 
B’s are known, then X;' will have a constant failure rate 4,°' and the procedures 


described in Section A can be used to obtain R(t) t(a) With T; replaced by 7," in 


equation (2.13). 
Suppose 8; is unknown and Xj7) , Xj (2), ---, Xi) are the ordered failure times 
under either type of truncated testing for component 7 in the system. Solutions 


B. and X, for B; and >; in the two equations given in equation (2.17) are the 


maximum likelihood estimates for B, and 4;. See Mann and others [Ref.1 pp 189-191]. 


These equations are used for both types of test truncation. If for component 1, 


iM 


testing is terminated on the f/" failure, then ¢,, = Xj) IN equation (2.17). The 
solution, B. , is a biased estimator for B;. Bain [Ref.5 pp 220] provides a table of 


constants B(n;) which depends on number of test items 7; such that e. = BB(n,) 


is a nearly unbiased estimator for 8. 
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If the testing for component ij is terminated at failure or at a given time ¢,; for 
each of the n; items on test, then ¢,, = t,; in equations (2.17a) and (2.17b). 

Now, let | 


ea with f= 


ij ij 2, uy Kk on (2.18) 
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In this thesis, the distribution of 7;, is approximated by the exponential distribution 


: : i * tee : : 
with failure rate ; =; and procedures similar to those in Section A are used 


to obtain the lower confidence limit on system reliability. Define 


fees (2.19) 


mete =) T, , t= 1, 2,..,k. Let i), = maxg, i. 


jal 


Note that this defines the index m. Define fF, by 





di T : 
aa — = ve meds eee (2.20) 
tore ai 


for both types of test truncation plans. 


Then an approximate 100(1-e)% upper confidence limit for 2, is given by 


mU(a) ~ —% . Gay 
2y° AT 
», tos 
where 
> f if test till f,"" failure 
F* = = ’ for all components + (2.22) 


14 > f if test till specified time 
"i? — for all components 


The corresponding approximate 100(1-a)% lower confidence limit R,(t),/,) for the 


reliability R,(t) of a series system is given by 


k : 

A * . B; 

R s(f)r (a) = expt = aaresps ore } wos (QO 
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when at least two components have at least one failure. Equations (2.11) and (2.12) 

also apply here when the total failures over all components is either zero or one. 
The accuracies of these approximate confidence interval procedures were 

evaluated by using computer simulations which are described in the next chapter. 


During this evaluation process, the degrees of freedom in the expressions X*, ,- 
in equation (2.21), X*,,, in equation (2.9) and  X*,5(;,.) in equation (2.7) were 


increased and decreased from the defined values of F and F given by these 
equations. The purpose of these modifications was to find more accurate lower 
confidence limit procedures. The specific increases and decreases are described in 
Chapter III. The results show that for some cases the procedures with modified 


degrees of freedom are more accurate. 
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lil. COMPUTER SIMULATION 


A. Test Plan 1: Testing n; Until f; Failures (RETP1) 

RETP1 is a program written in FORTRAN, on the Amdahl mainframe 
computer, which performs the computer simulation of the random failure times of 
the different types of components in the system. A documentation of this program 
and its associated subroutines is included in Appendix B. 

The program accepts input parameters via an input file IN1.DAT. For each 
replication, it generates the failure times for all the component items included in the 
test plan using a uniform random number generating subroutine LRNDPC. A quick 
evaluation of LRNDPC (see Appendix D) by plotting U(n+J) vs U(n) illustrates the 
uniformity of the routine. The program determines the total test time accumulated 
for each component in the system and computes the estimates of the key parameters 
and the consequent lower confidence limit for system reliability for that replication. 
The process is repeated 1000 times. When all replications are done, the routine 
EVAL processes the lower confidence limit estimates from all 1000 replications and 
determines the two measures of accuracy for the run, namely RSLOW and LEVEL. 

RSLOW is the 100(1-a@) percentile of the ordered set of lower confidence limits 
from the 1000 replications computed in a run. The true reliability of the system is 
RS. The closer RSLOW is to RS, the greater the accuracy of the procedure under 
evaluation in the run. If the procedure is exact, RSLOW will be equivalent to RS. 
To be conservative, RSLOW should always be lower than RS. 

LEVEL measures the proportion of 1000 lower confidence limits, from a run 
with 1000 replications, which are Jower than the true system reliability RS. The 
closer LEVEL is to the specified confidence level for the procedure, 1-a, the better 
the procedure. Values of LEVEL greater than (1-q@) reflect an under-estimation of 
RS which is conservative. Values of LEVEL lesser than 1-a signal an over- 


estimation of RS which may be undesirable. 


1a) 


Simulation runs are performed using RETP1 for all combinations of failure 


time distributions and levels of key input parameters listed below. 


(a) System. - 8 Exponential components in Series (Case 1) 
- 8 Weibull components in Series (Case 2) 
- 4 Exp and 4 Wei (Mixed) components in Series (Case 3) 


(b) True System Reliability (RS). 


- Hi (greater than 0.9) (Type A) 
- Lo (greater than 0.8) (Type B) 


(c) Level of Significance (a). 


- 0.1 
- 0.2 


(d) Degrees of Freedom for y? statistic (DF) as a function of the total number of 
failed test components (NFC) and total number of system components (NCOMP),. 


- DF = 2* NFC 

- DF = 2* (NFC + NCOMP ) 
- DF = 2 * (NFC-NCOMP ) 
- DF = 2 * NFC - NCOMP 


e) Test Plan for each component. 


- Test 5 until 5 failures 

- Test 15 until 15 failures 
- Test 15 until 11 failures 
- Test 15 until 7 failures 
- Test 15 until 3 failures 


For the 8 exponential components in Case 1, the mission time for each of the 
component is chosen to be 10 hrs. The program will accomodate different 
component mission times. The chosen values of the scale parameters, A; , were 
different depending on whether the system is highly reliable (Type A) or one with 
a lower reliability (Type B). The ratios between the largest and the smallest failure 
rate was chosen to be 8 and 4.5 respectively for Type A and Type B systems. 
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For the 8 Weibull components in Case 2, the mission time for each of the 
components was chosen to be 10 hrs. The chosen values of the scale parameters, X,, 
were different depending on whether the system is highly reliable (Type A) or one 
with a lower reliability (Type B). The ratio between the largest and smallest failure 
rate was chosen to be 8 for both system types. The shape parameter is chosen to be 
2 for all cases. The program will accomodate any value greater than zero for the 
shape parameter. 

A mixture of exponential and Weibull components with those parameters 
described in the last two paragraphs is chosen for the Type A and Type B systems 
of Case 3. 

Each simulation run of 1000 replication results in an output file OUT1.DAT. 
The raw output from all the RETP1 runs are summarized in tabular form and placed 
in Appendix E. Each table corresponds to a specific run case and system type 


combination. 


103. 


B. Test Plan 2 : Testing for a Specified Planned Test Time (RETP2) 

RETP2 is another program written in FORTRAN, on the Amdahl mainframe 
computer, which performs the computer simulation of the random failure times of 
the different types of components in the system. A documentation of this program 
and its associated subroutines is included in Appendix C. 

The structure of this program is quite similar to that of RETP1 described in 
Section A of this chapter. The program accepts input parameters via an input file 
IN2.DAT. For each replication, it generates the failure times for all the component 
items included in the test plan using LRNDPC. The program then determines the 
number of failed test components for each component in the system and computes 
the estimates of the key parameters and the consequent lower confidence limit for 
system reliability for that replication. The process is repeated 1000 times. When all 
replications are done, the routine EVAL processes the lower confidence limit 
estimates from all the 1000 replications and determines the two measures of accuracy 
for the run, namely RSLOW and LEVEL. The definitions of these two measures 
were discussed in the Section A. 

Simulation runs are performed using RETP2 for all combinations of failure 


time distributions and levels of key input parameters listed below. 


(a) System. = - 8 Exponential components in Series (Case 4) 

- 8 Weibull components in Series (Case 5) 

- 4 Exp and 4 Wei (Mixed) components in Series (Case 6) 
(b) True System Reliability (RS). 


- Hi (greater than 0.9) (Type A) 
- Lo (greater than 0.8) (Type B) 


(c) Level of Significance (qa). 


- 0.1 
- 0.2 
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(d) Degrees of Freedom for x’ statistic (DF) as a function of the total number of 
failed test components (NFC) and total number of system components (NCOMP). 


- DF = 2*(1+NFC) 
-DF = 13*2*(1+ NFC) 


(e) 10 values of K = 0.25, 0.5, 1, 2, 3, 4, 5, 10, 20, 30 


where K is a factor such that the expected number of failures for an exponential 
component during the specified total test time for that component is 0.6 times K. 
This K factor determines the bounds of the expected total number of failed test items, 
E[NFC]. The accuracy of the lower confidence limit procedures are highly correlated 
with E[NFC]. 

For the 8 exponential components in Case 4, the mission time for each of the 
component is chosen to be 5 hrs. The program can accomodate different component 
mission times. The chosen values of the scale parameters, \4,, were different 
depending on whether the resultant system is highly reliable (Type A) or one with 
a lower reliability (Type B). The failure rate was chosen to be 0.001 and 0.005 
failures/hr for all the components respectively. Total test time to be accumulated 
by each component is computed according to the following method. For each 
component 1 , t; represents the amount of operating time required to result in a 40% 
survival probability, that is, an expected failure of 0.6 component with an exponential 


failure time distribution. The computation for ¢; is as follows: 


R(t;) = exp(-At,) 


0.4 
] 
Gy In ( 0.4 ) 


e 
es 
tl 


T, , the total amount of test time to be accumulated for component! would be 
K times ¢; which will result in an expected number of 0.6 times K failed items for this 
component. E[NFC] will then be 8 times that number, since there are 8 such 


components in the system. 


Is 


For the 8 Weibull components in Case 5, the mission time for each of the 
component is chosen to be 15 hrs. The chosen values of the scale parameters, X. , 
were 0.005 failures/hr for the Type A system and 0.01 failures/hr for a Type B 
system. The shape parameter is chosen to be 2. For each Weibull component, a 
maximum of 20 test items are tested to failure. Estimation of E[NFC] and thus the 
total test time to be accumulated for each component is based on the exponential 
failure time model described in the earlier paragraphs. 

A mixture of exponential and Weibull components with those parameters 
described in the last two paragraphs is chosen for the Type A and Type B systems 
of Case 6. 

Each simulation run of 1000 replication results in an output file OUT2.DAT. 
The raw output from all the RETP2 runs are summarized in tabular form and placed 
in Appendix F. Each table corresponds to a specific run case and system type 


combination. 
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IV. RESULTS AND DISCUSSION 
The results of the simulation runs are summarized and discussed in this section. 
Tables 1A, 1B, 2A, 2B, 3A and 3B in Appendix E, and Tables 4A, 4B, 5A, 5B, 6A 
and 6B in Appendix F present the accuracy results in tabular form for all run cases 
that were simulated: a few of these tables appear in this section to facilitate 


discussion of the results. 


A. Test Plan 1: Testing m; Until f; Failures (RETP1) 

Table 1A displays the simulation results for Case 1 Type A. In this case, the 
system is comprised of 8 components in series. The failure time of each component 
has an exponential distribution. The failure rates of the 8 components range from 
0.0002 failures/hour to 0.0016 failures/hour. The mission time of the system is 10 
hours and the mission operating time of each component was also set to 10 hours. 
The component mission times do not need to be equal to the system mission time 
for the procedures evaluated in this thesis. This was discussed in Chapter II and is 
allowed for in all of the lower confidence limit equations. System reliability, RS, in 
Table 1A is 0.931 . Throughout this case, all of these parameters remain fixed. 

In simulation number 1 (S/N: 1) in Table 1A, five items for each of the 8 
components in the system are tested until they fail. Thus the number of failed 
components (NFC) is 40. One set of this 40 failure times is randomly generated for 
each simulation run. For this set of data, four 90% lower confidence limits and four 
80% lower confidence limits are computed. The four limits correspond to the values 
assigned to the degrees of freedom parameter F in the symbol x 2 2F which is a factor 


in the upper confidence limit equation for Ay.) .. The four different methods for 


computing this degree of freedom appear in the "Deg of Freedom" column. NCOMP 
denotes the number of components in the system. Thus for each simulation run, 


eight lower confidence limits are computed. After 1000 replications of this 


7 


Table 1A: 8 Exp in Series, RS = 0.931 (Hi) 
min » = 0.0002 f/hr, max X = 0.0016 f/hr, UT = 10 hrs 


Test Deg of eas Measures of Accuracy 







S/N Plan Freedom 








Test 5 until 
5 failed. 







NFC =40 






Test 15 until 
15 failed. 







NFC = 120 


3 | Test 15 until 2*NFC 0.927 0.955 


11 failed. 


NFC=88 
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Table 1A: 8 Exp in Series, RS = 0.931 (Hi) (Cont...) 
min A = 0.0002 f/hr, max X = 0.0016 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
RSLOW LEVEL 
2°NFC 
(112) 
0.923 0.931 
Pe os | em 
ee 















Test 15 until 
7 failed. 










NFC=56 











Test 15 until 
3 failed. 


NFC=24 


simulation are run, the 2 measures of accuracy RSLOW and LEVEL are computed. 
The lower confidence limit procedures are exact if RSLOW = RS in which case 
LEVEL = 1-a. 

Table 1A displays the accuracy results for 5 different sampling plans which are 
described in S/N: 1, 2, 3, 4 and S. 

A comparison of the four values of RSLOW for each of these five sampling 


plans reveal that the lower confidence limit procedure with degrees of freedom equal 
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to 2*NFC-NCOMP is the most accurate lower confidence limit procedure. In S/N 
1, for example, the RSLOW value of 0.927 is the largest such value below the RS 
value of 0.931. Values of RSLOW above RS are optimistic and not as desirable as 
values of RSLOW which are equi-distant below RS. 

The values of RSLOW and LEVEL are based on 1000 replications and their 
accuracy merit should roughly be measured to the nearest one hundredth. That is 
we should round 0.927 to 0.93 and compare it with RS = 0.93. It is evident that the 
lower confidence limit procedure with degrees of freedom equal to 2*NFC-NCOMP 
is very accurate for all cases simulated. 

Table 2A displays the accuracy results of Case 2 for Type A systems. In this 
case, the 8 components connected in series have the shape parameter 8 = 2 and the 
scale parameters i, varying between 0.001 and 0.008 failures/hr. Mission time is 10 
hours and each component has this same mission or utilization time (UT). 
Inspection of Table 2A reveals the following: 

(1) More than 5 items of each component should be tested until failure for 
any of these procedures to be reasonably accurate. 

(2) If 15 items of each component are tested until all fail, then these 
procedures will be reasonably accurate when the degrees of freedom is either 
2*NFC-NCOMP or 2*(NFC-NCOMP),. 

(3) The procedures are reasonably accurate for 80% confidence level when 
the truncation is not below 7 out of 15 items. 

(4) The procedures are slightly conservative at the 90% confidence level when 
the truncation is 7 out of 15 items or 11 out of 15 items. 

There are numerous ways to modify these lower confidence limit procedures 
to effect improvements in their accuracy. One avenue is to modify the estimate for 
the shape parameter, 8. Some very recent work in the literature provides a method 
for estimating 8 that differs greatly from the maximum likelihood estimator (MLE) 


and does not require computer iteration. 
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Table 2A : 8 Wei in Series, RS = 0.980 (Hi) 
min X = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 


Test Deg of Measures of Accuracy 
Plan Freedom 


RSW | SE VeEL 
a [avec | on | ose [oom | 
Toa | ovo | one 













= 





Test 5 until 
5 failed. 


a 






NFC=40 





0.937 0.994 


[ay ree BE re 


0.933 0.981 
2"NFC 0978 
(240) 
0.2 0.974 0.913 


0.977 0.931 
iam Sam 


3 Test 15 until 2*NFC 0.982 0.876 
11 failed. 
0971 


° 090 
pesos | ow 


0.2 0.979 0.819 


Test 15 until 
15 failed. 






NFC = 120 





NFC=88 


= 
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Table Ze 8 Wei in Series, RS = 0.980 (Hi) (Cont...) 
min X = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 


Test Deg of Measures of Accuracy 
S/N Plan Freedom 
RSLOW LEVEL 


Far ae ee = 






7 failed. (112) 


NFC=56 


Test 15 until 
3 failed. 


NFC=24 


Table 3A displays the results of Case 3 for Type A systems. In this case 8 
components are connected in series. Four of them have failure times with 
exponential distributions and the remaining four have failure times with Weibull 


distributions each with shape parameter, 8, equal to 2. 


IL 


Table 3A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) 
min A = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs 





Test Deg of 
S/N Plan Freedom 
1 


Test 5 until 2*NFC 0.979 0.942 
5 failed. 80 
Se 


NFC=40 





Test 15 until 
15 failed. 


NFC = 120 





3 | Test 15 until 2*NFC 0.986 
11 failed. (176) 
fai im 


NFC=88& 





| | 2*(NFC- 0.982 0.702 
| NCOMP) — 
160 0.2 0.982 0.591 


Ma 


Table 3A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) (Cont...) 
min = 0.002 f/hr, max X = 0.008 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
RSLOW LEVEL 


Teas | =r 
ie eer 


NFC=56 = 0.982 0.952 


0 978 


Test 15 until 
3 failed. 


NFC=24 2e(NEC + 





Inspection of Table 3A reveals the following: 

(1) The two procedures corresponding to degrees of freedom equal to 2*NFC 
and 2*NFC-NCOMP are reasonably accurate for all 5 simulation cases. 

(2) The procedures appears to be nearly equally accurate for both 80% and 
90% confidence levels. 
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B. Test Plan 2: Testing for a Specified Planned Test Time (RETP2) 

In the simulations for test plan 2, components were tested until failure or until 
some planned test time scenario. Failed items were not replaced. Components 
whose failure times had exponential distributions were tested until a pre-determined 
total test time was accumulated for their type of component. Components whose 
failure times had Weibull distributions were tested until failure or a pre-determined 
planned test time for that test item. The latter truncation plan is needed for Weibull- 
type items in order to use the maximum likelihood estimates [as in equation (2.17)] 


to solve for £ . 


Inspection of Table 4A reveals that the lower confidence limit procedure for 
degrees of freedom equal to 2*(1+NFC) is quite accurate when enough testing is 
done to make the expected number of failures, E[NFC], greater than or equal to 4.8. 
This testing constraint is well within the domain of constraints set on testing in 
development programs for major systems within the Department of Defense. 

Examination of Table 5A reveals that the lower confidence limit procedure for 
degrees of freedom equal to 2*(1+ NFC) is moderately accurate when enough testing 
is done to make E[NFC] greater than or equal to 9. The accuracy diminishes slightly 
as E[NFC] increases. This could be corrected by decreasing the degrees of freedom 
slightly to make RSLOW slightly larger. 

The results displayed in Table 6A show that the lower confidence limit 
procedure for degrees of freedom equal to 2*(1+NFC) is quite accurate when 


enough testing is done to make E[NFC] greater than or equal to 9.6. 


2D 





















Table 4A : 8 Exp in Series, RS = 0.961 (Hi) 
= 0.001 f/hr, UT = 5 hrs 
Degrees See 
S/N of (TT) 
ae 


1 | 2a+NFC) 
05 / 24 
= 

| 10/48 

= 

20 | 96 

(1800) | ee 

~ 
so/2 | o1 | ose | ass 
ae 

100 / 48 

a 

210 | 96 

= Fea | see ea 


30.0 / 144 0.906 
(27000) 
0.2 0.961 0.804 
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Table 5A : 8 Wei in Series, RS = 0.956 (Hi) (*) 
= 0.005 f/hr, UT = 15 hrs 























Degrees K / E[NFC] Measures of Accuracy 
of (TT) 
Freedom {| RSLOW LEVEL 


— a I - a _ SS 


1 | 2aeNFo) | 025/12 i | 
4 
20 / 96 

ad 

5.0 / 24 
em 

20.0 / 96 


30.0 / 144 0.1 0.952 0.940 


(5400) 
02 | 0.946 0.928 












(*) 20 test items for each Weibull component. 
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Table 6A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) 
(exp) = 0.001 f/hr, UT(exp) = S hrs 
(wei) = 0.005 f/hr, UT(wei) = 15 hrs 


7 
eee RSLOW LEVEL 
(225) 
(5) ics ec 
0.957 0.803 
(2700) 
S40 ass | ost 
(3600) 
(720) 0.957 0.839 
| 


(27000) 
5400 0.956 ae 862 
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The accuracy results of simulations performed in this thesis cannot be extended 
to systems that differ significantly from those simulated here. However, it can be 
said that the procedures which are accurate for series systems are also usually 
accurate for 1-out-of-k parallel systems because system reliability, R,, can be written 


in terms of component reliabilities, R;, as 


k 


R,-1-[] (1 - R,) 


i-l 
k 
Thus, upper confidence limits on [[ (1 - R;) will yield a lower confidence 


i=] 


limit on R, The accuracy of the upper confidence interval procedures for 
k 


[] (1 - 8,) ,aseries-type problem, should be nearly the same as those obtained 
i=1 


in this thesis because equations like (2.9) would be replaced with equations for the 


lower confidence limit 4,,)(,) ON A, and would look like 


2 
. X 1-a,2F 


m,L (a) 7 k 
Dyan, 


i=] 


. (4.1) 


If the degrees of freedom parameter F is large, as it is in the cases simulated in this 
thesis, the associated Chi-Square distribution with 2F degrees of freedom is nearly 
symmetric about its mean so the lower tail has nearly the same shape as the upper 


tail. This characteristic should yield very similar accuracies for i, l(a) aS was 
obtained for ares in this thesis. These accuracy comparisons translate directly 


to similar comparisons about the accuracies of the associated lower confidence limit 


procedures for the system reliability of series and parallel systems. 


ay 


Table 7 : Summary of Procedure Accuracy by Simulation Category (RETP1) 


S/N Key Parameters and Observation and Discussion 
their Levels based on RSLOW and LEVEL 


= — a a a a a 







1 Run Cases = For all exponential systems, accurate | 
» All Exponential (1) procedures were developed for component | 
s All Weibull (2) Sample sizes > 5. 

= Mixed (3) = For all Weibull systems, component 
sample sizes should be > 15 with 
truncation at r > 7 failures. 

































System Reliabilities 
= Hi (> 0.9) (Type A) 
» Lo (> 0.8) (Type B) 


Accurate procedures were developed for 
both cases if sample sizes are adequate. 


Levels of Significance 
na = 0.1 
aa = 0.2 


Accuracy varied slightly depending on 
system type and test plan, but accurate 
procedures exist for both levels. 






















Degrees of Freedom 
s DF = 2*NFC 

» DF = 2*(NFC+NCOMP) 
a» DF = 2*(NFC-NCOMP) 
a DF = 2*NFC-NCOMP 


Greatest accuracy as follows: 
s All exponential : 2* NFC-NCOMP 
s All Weibull : 2* NFC-NCOMP 

or 2*(NFC-NCOMP) 












Test Plan 
s Test 5 until 5 failures 

s Test 15 until 15 failures 
s Test 15 until 11 failures 
s Test 15 until 7 failures 
» Test 15 until 3 failures 


Accurate procedures existed for all run 
cases for all system types except for the all 
Weibull system where number of test 
items n > 15 and number of failures r 
should be > 7. 
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Table 8 : Summary of Procedure Accuracy by Simulation Category (RETP2) 


» All Exponential (4) 2*(1+NFC) when enough testing was 
» All Weibull (5) done to make the expected number of 
» Mixed (6) failed components E[NFC] > 9. 


System Reliabilities Same as S/N(1). 
a Hi (> 0.9) (Type A) 
a Lo (> 0.8) (Type B) 


Levels of Significance Same as S/N(1). 
aa = 0.1 
aa = 0.2 


Degrees of Freedom DF = 2*(1+NFC) was the most accurate 
» DF = 2*(1+NFC) procedure. 
a DF = 1.3*2*(1+NFC) 


Test Plan K should be chosen so that E[NFC] > 9. 
» K factors of 

0.25.0, 2,03, 

4,5, 10, 20, 30 





Tables 7 and 8 provide cursory summaries of some constraints needed to assure 
the existence of one or more accurate lower confidence limit procedures among the 
procedures that were evaluated. The simulation scenarios are divided into five 


categories for this summarization. 
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V. APPLICATION EXAMPLES 
Based on the procedures evaluated by the RETP1 and RETP2 runs, four 
different test plans and failure time data were constructed to illustrate the use of the 
procedures in providing a lower 100(1-a) % confidence limit for the system reliability 


of a series system with different types of components. 


CASE 1: 8 Exponential Components in Series 
=SSSSes TEST PLAN 1 - Test 15 until 7 fails for each component 


I. Raw Data 
Comp Ordered Failure Times (h) 
t sae) T(2) T(3) T(4) T(5) T(6) T(7) 
i 300.0 400.0 500.0 600.0 700.0 800.0 900.0 
2 350.0 450.0 550.0 650.0 750.0 850.0 950.0 
S 400.0 500.0 600.0 700.0 800.0 900.0 1000.0 
4 450.0 550.0 650.0 750.0 850.0 950.0 1050.0 
5 500.0 600.0 700.0 800.0 900.0 1000.0 1100.0 
6 550.0 650.0 750.0 850.0 950.0 1050.0 1150.0 
7 600.0 700.0 800.0 900.0 1000.0 1100.0 1200.0 
8 650.0 750.0 850.0 950.0 105059 1150.0 1250.0 
II. Data Summary 
Comp ER(i)* 


1 5.0 TS 7 11400.0 0.00053 1.00000 11400.0 
2 5.0 IG 7 12150.0 0.00049 0.93827 11400.0 
3 320 lis} 7 129800.0 0.00047 0.88372 11400.0 
4 aah) Ls 7 13650.0 0.00044 0.83516 11400.0 
5 5.0 rs 7 14400.0 0.00042 0.79167 11400.0 
6 5.0 1S 7 15150.0 0.00040 0.75248 11400.0 
u 3.0 ales 7 15900.0 0.00038 0.71698 11400.0 
8 5.0 15 7 16650.0 0.00036 0.68468 11400.0 
III. Estimation Procedure for RSLOW 
Parameter df Value 
ALPHA 0.1 
NFC 56 
CHISQD lal: 131.56 (df = 2 * NFC , CHISQD from tables) 
LMU 0.00072 
RSLOW 0.97647 
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CASE 2 : 8 Weibull Components in Series 
<7" TEST PLAN 1 - Test 15 until 7 fails for each component 


I. Raw Data 
Comp Ordered Failure Times (h) 
sh T(1) T(2) T(3) T(4) ECS) T(6) T(7) 
al 10.0 20.0 20.0 40.0 50.0 60.0 70.0 
2 20.0 30.0 40.0 50.0 60.0 70.0 80.0 
3 30.0 40.0 50.0 60.0 70.0 80.0 $0.0 
4 40.0 50.0 60.0 70.0 80.0 90.0 100.0 
5 207.0 60.0 70.0 80.0 90.0 100.0 110.0 
6 60.0 70.0 80.0 $0.0 100.0 £1020 120.0 
7 70.0 80.0 $0.0 100.0 110.0 120.0 130.6 
8 80.0 90.0 100.0 110, 0 120.0 130.0 140.0 
II. Data Summary 
Comp ER(i)* 


i UTC.) NC(i) NF (i) TT(i)  ELM(i) ER(i) ml CL) 


wee eee eee er eer eae ewe SO ew were ee 2B ee 2 S28 @e FF Fe FT F@Beesw swe’ ee FZ Fee eres ZT SO SF S88 88 SF SOF FT 8 Te sBe se see =e = — 








1 520 15 7 5.3E+04 1.32E-04 1.00000 5.3E+04 
2 5.6 15 7 7.2E+04 9.7SE-05 0.74406 5.3E+04 
3 5.0 15 7 9.3Et+04 7.54E-05 0.57328 5.3E+04 
4 570 15 7 1.2E+05 5.98E-05 0.45431 5.3E+04 
5 5.0 a5 7 1.4E+05 4.85E-05 0.36842 5.3E+04 
6 5.0 5 7 1.7E+05 4.01E-05 0.30452 5.3E+04 
7) 5.0 15 7 2.1E+05 3.37E-05 0.25577 5.3E+04 
8 5.0 15 7 2.4E+05 2.87E-05 0.21777 5.3E+04 
III. Estimation Procedure for RSLOW 
Parameter df Value 
BETA 2.0 
ALPHA 0.1 
NFC 56 
CHISQD 112 131.56 (df = 2 * NFC , CHISQD from tables) 
LMU 1.55E-04 
RSLOW 0.98497 
IV. Workarea 
Comp Ordered Failure Times (h) raised to the power of BETA 


= T’(1) T’ (2) T' (3) T’ (4) C3) T*(G) a IC?) 


ae —-S Ss = Se Sew ee ew KH ww ww ww www ww www ww ww ww www Ow OM Ow Ow Mw we SE SE SE SE SOO SO Se ST Ow SE SE Bw Ow ww ew www ew ew ewe Se 
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CASE 3 : 8 Exponential Components in Series 
==<==-==- TEST PLAN 2 - Test until TT(i) for each component 


I. Raw Data 


No data except noting the number of failures for each component (NF(i)) 


II. Data Summary 


Comp ER(i)*  ER(i)* 
i UT (i) NF (i) TT (i) ELM(i ) ER (i) TT(i) UT (i) 
ak a0 6 2000.0 0.00300 0.85714 1714.3 4.3 
2 5.0 6 2000.0 0.00300 0.85714 1714.3 4.3 
3 5.0 5 2000.0 0.00250 0.71429 1428.6 3.6 
4 5.0 7 2000.0 0.00350 1.00000 2000.0 5.0 
5 5.0 5 2000.0 0.00250 0.71429 1428.6 3.6 
6 5.0 5 2000.0 0.00250 0.71429 1428.6 3.6 
7 50 4 2000.0 0.00200 0.57143 1142.9 2.9 
8 5.0 7 2000.0 0.00350 1.00000 2000.0 5.0 
III. Estimation Procedure for RSLOW 
Parameter ob Value 
ALPHA 0 jee 
NFC 45 
CHISQD 92 109.76 (df = 2 * ¢( 1+ NFC ) , CHISQD from tables) 
LMU 0.00427 
RSLOW 0.87180 
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CASE 4 ; 8 Weibull Components in Series 
a — TEST PLAN 2 - Test 20 until TT(i) for each component 


I. Raw Data 
Comp Ordered Failure Times (h) 
i (4°) TCT) T(2) T(3) T(4) TG) T(6) 
1 60.0 10.0 20.0 30.0 40.0 50.0 60.0 
2 60.0 15.0 25.0 35.0 45.0 55.0 
3 60.0 yAt a8) 30.0 40.0 50.0 60.0 
4 60.0 25.0 35.0 45.0 55.0 
> 60.0 30.0 40.0 50.0 60.0 
6 60.0 a5 0 45.0 55...0 
7 60.0 40.0 50.0 60.0 
8 60.0 45.0 220 
II. Data Summary 
Comp ER(i)* 





1 5.0 20 6 6.0E+04 1.01E-04 1.00000 6.0E+04 
2 5.0 20 5 6.1E+04 8.18E-05 0.81118 5.0E+04 
3 5.0 20 5 6.3E+04 7.94E-05 0.78704 5.0E+04 
4 50 20 4 6.5E+04 6.20E-05 0.61499 4.0E+04 
5 5.0 20 4 6.6E+04 6.04E-05 0.59919 4.0E+04 
6 5.0 20 3 6.7E+04 4.45E-05 0.44090 3.0E+04 
7 5.0 20 3 6.9E+04 4.35E-05 0.43179 3.0E+04 
8 5.0 20 2 7.0E+04 2.86E-05 0.28394 2.0E+04 
III. Estimation Procedure for RSLOW 
Parameter df Value 
BETA 2.0 
ALPHA (3 
NFC 32 
CHISQD 66 81.09 (df = 2 * ( 1+ NFC ) , CHISQD from tables) 
LMU 1.28E-04 
RSLOW 0.98425 
IV Workarea 
Comp Ordered Failure Times (h) raised to the power of BETA 


i TT’ (1) ae T) Tez) (3) T’ (4) TCS) T’ (6) 


1 3.6E+03 1.0E+0O2 4.0E+O02 9.0E+O2 1.6E+03 2.5E+03 3.6E+03 
2 3.6E+03 2.3E+OZ 6.3E+O0Z2 1.2E+03 2.0E+03 3.0E+03 

3 3.6E+03 4.0E+O2 9.0E+O2 1.6E+03 2.5E+03 3.6E+03 

4 3.6E+03 6.3E+t0O2 1.2E+03 2Z.0E+03 3.0E+03 

> 3.6E+03 9.0E+O2 1.6E+03 2.5E+03 3.6E+03 

5 s.0nt0s 1.ZEt03) Z.0EtTOS 3 .0ET0S 

7 3.6E+03 1.6E+03 2.5E+03 3.6E+03 

8 3.6E+03 2.0E+03 3.0E+03 


35 


VI. CONCLUSION 
Some of the lower confidence limit procedures developed and evaluated in this 
thesis are reasonably accurate for the series systems simulated for test plans with 
sample sizes and truncation scenarios that are usually experienced in DoD aquisition 
programs. ‘The accuracy of these methods can be varied by modifying the degrees 


of freedom parameter, F, in the x’, ; term in equation (2.9), namely: 


2 

Amite) = 2. 

2d rT. 

The computer program can be modified with modest effort to accomodate specific 
complex coherent systems so long as the components have failure time distributions 
that are exponential or Weibull. This means that the computer program provided in 
this thesis can be used to develop a reasonably accurate lower confidence limit for 
the system reliability of a specific complex quasi-coherent system with independent 
components. This can be done by choosing the failure distribution and associated 
parameters of the components, the corresponding test plan parameters and the 
desired level of confidence. The simulation can then be run for this set of 
parameters for various equations for the degrees of freedom parameter, F, to 
determine an equation for F that yields a lower confidence limit with a satisfactory 
degree of accuracy. 

When testing is truncated on the number of failures, a reasonably accurate 
procedure existed for all 3 cases of systems (all exponential, all Weibull and mixed) 
that were simulated when (a) the sample size of the components was 10 or larger, 
and, (b) the ratio of the number of failures to the sample size was at least 0.5. When 
testing was truncated by planned test time, reasonably accurate procedures were 


found for cases where the expected number of failures was at least 7. 
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VII. RECOMMENDATIONS 
The computer program developed in this thesis facilitates the development of 
lower confidence limit procedures for explicit quasi-coherent systems. Systems other 
than series systems with large numbers of components (eg. 30) should be simulated 
to test the versatility of the general lower confidence limit methods used here. 
Modified estimates for the 8 (shape) parameter in the Weibull failure time 
distribution and the parameter r = i,/X,, should be explored in an attempt to find 


more accurate procedures. 


ay) 


APPENDIX A: Derivation of Formula Used 


Suppose T has an exponential distribution with failure rate 4. Suppose Ty, Ty» 


T (3) » Tq) are the first r ordered statistics in a random sample of size n from this 


exponential distribution. Let S be defined by 


S = » Ty + (n-r)Ty 


It is well known that 2XS has a Chi-Square distribution with 2r degrees of freedom 
(See Ref.3 p 488). The maximum likelihood estimator for ) is given by 


ae 
5 


We seek an unbiased estimator for A . Suppose X has a Chi-Square distribution 


with 2r degrees of freedom, then 


. = 1 eal _x 
IOS) PT) exp( ~~ 


and the integral of this function from zero to infinity equals 1 (See Ref.3). 


E(=) = { <x Lexp(-2) dx 


J 2’ T(r) 
_T(r-1)p 1 (1) eee 
2T(r) Job C=) mK 2? 
| 
; Fe) 
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Then 


1 ae rd 
XS (r-1) 





E(3) = E(=) = r2E( 


Therefore 


E(—=3) : E() =) 


BY 


APPENDIX B: Users’ Guide for RETP1 


Reliability Estimation Test Plan 1 (RETP1). 
by YEE, Kah-Chee 


July 91 





1. Brief Description. 


RETP1 is a computer program written in FORTRAN that runs on the Amdahl 
mainframe at NPGS. It allows the user to simulate exponential and Weibull failure 
times of component items being tested to evaluate the accuracy of a confidence limit 
estimation procedure based on Type II data censoring (that is, testing n, items of 
component i until f; of them failed). | 


2. Program Input. (IN1.DATA) 


The input of the program are specified to the program via an input file called 
IN1.DAT. A sample input file is shown below. 


This file contains the inputs required by the RETP1 model. 
Update only the numerical values between dotted lines as appropriate. 


Do not delete any of the comment lines. (INL.DAT) 14 May 91 
een mens - A et 
C Value Type Units Description Variable 
en enn ny ee a a ee” 
16807 .0 Real = initial random seed ISEED 
8 Int - total # of components in system NCOMP 
4 Int - # of EXPonential components NEXP 
4 Int - ## of WEIbull components NWEI 
0 Int - # of GEOmetric components NGEO 
GO.08 Real ° tolerance for MLE TOL 
0.20 REAL - DESIRED SIGNIFICANCE LEVEL ALPHA 
1000 Int - # of replications desired NREP 
3 Int - test case number TCN 
1 = all EXP 
2 = all WEI 


3 = EXP + WEI 
4 = EXP + WEI + GEO 
8 Int - number of cut sets NCS 


C TEST PLAN : Testing NC(I) items of component i 
© until NF(I) of them fails. (Use REAL numbers ONLY!!!) 
nena ee Oe ee ee ee e-- 
C Comp Comp Comp Parameters Util Test Plan Inputs 
C Number Type Scale Shape Time/Cycle # Comp # Failed 
c WT TY(I) PARM1(I) PARM2(1) UT(I) UC(CI) NC(I) NF(I1) 
6 Int Int Real Real (hrs) Int Int Int 
Inner ete rere less eo aia 2 else oo eS Se eS See et ee ee ee ee 
1.0 bd 0.0020 1.0 10.0 15.0 3.0 
20 0 0.0040 1.0 10.0 on 3.0 
BO 1.0 0.0060 ro 10.0 15.0 3.0 
4.0 1.0 0.0080 8, 10.0 15.0 3.0 
5.0 ZO 0.0020 2.0 10.0 15.0 320) 
6.0 2.0 0.0040 220 10.0 1330 3.0 
7.0 2.0 0.0060 2.0 10.0 LaaO 3.0 
8.0 2.0 0.0080 jam, 10.0 15.0 3.0 
ble. oe See 
C Note : TY(I)=1 EXPONENTIAL P(surv) = exp(-PARM1)*T) 
C TY(I)=2 WEIBULL P(surv) = exp(-(PARM1*T)**PARM2) 
C TY(I)=3 GEOMETRIC P(surv) = PARM1**N 
Cae esee coco cd so GES Sa = ae 
C SYSTEM CONFIGURATION : Identification of CUT SETs 
C - min groups of components that have to fail 
C for the system to fail. 
sa eesccgees OOO Be I 6 ee 
C Cut Set ## in Set List of Components in Cutset 
C J COMP(J,1) COMP(J,2) ... up to COMP(J,1) components 
Lo a Se ge S SRS S S666 66 eee ee eee ee ee 
ih 1 ao 0 80-0 Oe Omer 0 
2 a Zo OSROeeO TO” 0. 0720 0 
3 ir 3 0 0 0 0 0 0 0 0 
4 iL 4 00 00 0 0 0 0 
5 Ht Spl al |) eas 0 @ a © mes © Tm ©. 
6 I Guwe OT OO. 0 OF 0 0 
i 1 7 0 0 0 0 0 0 0 0 
8 i Ba te Oe Oar OO: FONO® 0 
IEEE eS ee ns coe on co ce eee 6 6 6 3 Se + - eee 


3. Program Flow and Logic. (NAME1.DEF, PARM1.DEF and RETP1.FOR) 


Input parameters are first read in by the program by calling the INPUT 
subroutine. The program then evoke the SIM subroutine which generates the 
random failure times and compute the key statistics required in the procedure. The 
next subroutine EVAL determines the measures of accuracy for run. REPORT is 
the subroutine which generates the output file for the run OUT1.DAT. 


The variables in the program RETP1.FOR are described in the file 
NAME1.DEF as listed below. 
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This file contains the declaration for input and output variables 
used in the the RETP1 model. (NAME1.DEF) 14 May 91 


ISEED = initial random seed selected. 

SEED = current random seed. 

RS = true overall series system reliability. 

ALPHA = level of significance desired. 

NREP = number of replications desired for the simulation. 
TPN = test plan number (1). 

TCN = test case number (1, 2, 3 or 4). ; 

NCOMP = total number of components in the system. 

NEXP = number of components with EXP failure times. 
NWEI = number of components with WEI failure times. 
NGEO = number of components with GEO failure times. 
TOL = desired tolerance for MLE of WEI shape parameter. 


Distribution: EXPonential WEIbull GEOmetric 

TY(I) = type: 1 . 3 
PARM(1,1): Scale(l1/hr) Scale(l/hr) Prob 
PARM(2,1): Shape - 


UT(I) = utilization time (hrs) for component i (EXP and WEI). 

UC(I) = utilization cycles for component i (GEO only). 

NC(I) = number of test samples (sample size) for cmponent i. 

NF(I) = desired number of failures in test for component i. 

NCS = number of cut-sets for the system. 

COMP(J,K) = kth parameter of cut-set j (first being the no. of 
components belonging to the cut-set). 


Assumed Variables. 

MAXCOMP = maximum number of components allowed in the system. 
MAXREP = maximum number of replications permitted. 

MAXCUT maximum number of cut-sets. 


Program and Output Variables. 
RS = true overall system reliability. 
TT(I) = total accumulated failure time (hr) for component i 
(EXP and WEI only). 
TC(I) = total accumulated cycles to failure (incl. failure cycle) 
for component i (GEO only). 
EBETA(I) = estimate for shape parameter of component i (if Weibull). 
REL1(J) = actual reliability for cut-set j. 
REL2(J) = computed reliability for cut-set j for current replication. 
ELM(I) = estimated component failure rate (l/hrs) for component i. 
ELMAX(M) = max estimated component failure rate for rep m (l/hrs). 
ER(I) = ratio of estimated failure rate to LMAX. 
NFC(M) = total number of failed test components. 
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LMU(M) = upper confidence bound for failure rate (l/hrs), 
RSL(M) = lower confidence limit estimated for system reliability 
for the mth replication. 
ORSL(M) = ordered RSL(M) (ascending). 
RSLOW = (1-ALPHA)x100 percentile of set of RSL(M). 
LEVEL = achieved confidence level, ie. proportion of RSL(M) that 
are lesser than RS (conservative estimate). 


FPO TYAS L Wise) cs 5 3 SS a ee 
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Together with the main program in RETP1.FOR are the other subroutines 
needed in the simulation. The declaration of variables is done in the file 
PARM1.DEF. Relevant descriptions are included as comment lines in the source 
code to help explain the program segments. A listing of PARM1.DEF and 
RETP1.FOR is given below. 


C This file contains the declaration for input and output variables 
C used in the the RETP1 model. (PARM1.DEF) 14 May 91 


INTEGER MAXCOMP, MAXREP 

PARAMETER( MAXCOMP = 100 , MAXREP = 1000 , MAXCUT = 20 ) 

REAL*8 ISEED, SEED 

INTEGER NREP, TCN, NCOMP, NEXP, NWEI, NGEO, NCS, 
NC(MAXCOMP), NF(MAXCOMP), TY(MAXCOMP), NFC(MAXREP), 
UC(MAXCOMP), TC(MAXCOMP), COMP(MAXCUT,MAXCOMP) 

REAL*8 RS, ALPHA, UT(MAXCOMP), TT(MAXCOMP), 


%~ + 


* PARM(2,MAXCOMP), ELM(MAXCOMP), ER(MAXCOMP), 
* LMU(MAXREP), RSL(MAXREP), ORSL(MAXREP) , 
* ELMAX(MAXREP), RSLOW, LEVEL, TOL, EBETA(MAXCOMP) , 
* REL1(MAXCUT) , REL2(MAXCUT) 
C 
COMMON/BLOCK1/ISEED, SEED, NREP, TCN, NCOMP, NC, NF, NEXP, NWEI, 
* NGEO, NCS, TY, NFC, UC, TC, COMP 
COMMON/BLOCK2/RS, ALPHA, UT, TT, PARM, ELM, ER, LMU, 
* RSL, ORSL, ELMAX, RSLOW, LEVEL, TOL, EBETA, 
* REL1, REL2 
c 
ce (S)NT QI JOM RIT DTS e oc GE SE a nn el eg 
C 
CG ey HENS See Se oe ce eee el 


C This file contains the main program and the subroutines 
C for the Reliability Estimation Test Plan 1 (RETP1) model. 
C (RETP1.FOR) - runs on a IBM PC Compatible. 

C 

C MAINFRAME VERSION... 

C 

C Test Plan 1 : Testing NC(I) items for component i 
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oe Until NF(I) of them fails. 


by Yee Kah-Chee SMC 2802. 
14 May 91. 


PROGRAM RETP1 
Include the declaration files. 


INCLUDE ‘NAME] DEF’ 
INCLUDE ‘PARM] DEF’ 


Read in input data. 
CALL INPUT 
Activate simulation. 
CALL SIM 
Process and evaluate output data. 
CALL EVAL 
Generate simulation report. 
CALL REPORT 


STOP 
END 


SUBROUTINE INPUT 
This subroutine reads in the inputs for the RETP1 model. 
Include the declaration file. 

INCLUDE ‘PARM]1 DEF’ 


INTEGER I, J, K, DUM2(11) 
REAL*8 DUM1(7) 


Read data from ‘IN1.DAT’ designated as logic unit l. 


OPEN(UNIT=1) 
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10 


20 


30 


50 


READ(1, 10) 
FORMAT (1X, /////) 
READ(1,*) ISEED 
READ(1,*) NCOMP 
READ(1,*) NEXP 
READ(1,*) NWEI 
READ(1,*) NGEO 
READ(1,*) TOL 
READ(1,*) ALPHA 
READ(1,*) NREP 
READ(1,*) TCN 
READ(1, 20) 
FORMAT(1X,///) 
READ(1,*) NCS 


READ (1,30) 
FORMAT (1X, ////////) 


DO 50 K = 1, NCOMP 


READ(1,*) DUM1 


I = NINT(DUM1(1)) 
TY(I) = NINT(DUM1(2)) 


IF (TY(I).EQ.1) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 
UT(I) = DUM1(5) 

NC(I) = NINT(DUM1(6)) 
NF(I) = NINT(DUM1(7)) 
EBETA(I) = 0 

ELSEIF (TY(1).EQ.2) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 
UT(I) = DUM1(5) 

NC(I) = NINT(DUM1(6)) 
NF(I) = NINT(DUM1(7)) 
ELSEIF (TY(1).EQ.3) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 


UC(I) = NINT(DUM1(5)) 
UT(I) = DUM1(5) 
NC(I) = NINT(DUM1(6)) 
NF(I) = NINT(DUM1(7)) 
ENDIF 
CONTINUE 
READ (1,60) 


60 FORMAT(1X,///////////) 


DO 80 I = 1, NCS 
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READ(1,*) DUM2 
J = DUM2(1) 
COMP(J,1) = DUM2(2) 
DO 70 K = 1, COMP(J,1) 
COMP(J,K+1) = DUM2(K+2) 
70 CONTINUE 
80 CONTINUE 
CLOSE(UNIT=1) 
RETURN 
END 


SUBROUTINE SIM 


This subroutine simulates NREP possible outcomes of the test plan 
desired in order to obtain the raw estimates of LMU(M) and RSL(M) 
for each of the replication. 


Include the declaration file 
and declare local variables. 


INCLUDE ‘'PARM] DEF’ 

INTEGER I, J, Kyu, 1SUM)) REx 
REAL*8 UNI 

REAL*8 SUM, PROD, 
* FT(MAXCOMP) , OFT (MAXCOMP) 


SEED = ISEED 
Compute overall true system reliability RS. 


RS = 1.0 
DO 30 J = 1, NCS 
PROD = 1.0 
DO 20 I = 1, COMP(J,1) 
K = COMP(J,I+1) 
PROD = PROD*( 1 - SURV(TY(K) ,PARM(1,K),PARM(2,K),UT(K)) ) 
20 CONTINUE 
REL1(J) = 1.0 - PROD 
RS = RS * REL1(J) 
30 CONTINUE 


Start of Simulation. 
(Initialize replication counter M). 


M=1 
DO WHILE (M.LE.NREP) 


PRINT 35, M 
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35 FORMAT(1X,’Replication ' ,I4) 


Test Plan : Sample and determine unknown TT(1) 
eee eee eee with known NC(I) until NF(I) fails. 


Generate NC(I) failure times, put them in ascending order 
with the smallest failure time on the top of the list. 


DO 70 I = 1, NCOMP 


DO 40 K = 1, NCCI) 
CALL LRNDPC(SEED,UNI,1) 


IF(TY(I).EQ.1) THEN 
FT(K) = -LOG(UNI) /PARM(1,1) 
ELSEIF(TY(1I).EQ.2) THEN 
FT(K) = (1.0/PARM(1,1))*(-LOG(UNI) )**(1.0/PARM(2,I1)) 
ELSEIF(TY(I).EQ.3) THEN 
FT(K) = 1.0 
DO WHILE (UNI.LT.PARM(1,1)) 
FT(K) = FI(K) + 1.0 
CALL LRNDPC(SEED,UNI,1) 
ENDDO 
ENDIF 
40 CONTINUE 


Bubble Sort the failure times in ascending order. 
CALL BUBBLE(NC(I) ,FT,OFT) 


Compute the total time accumulated in the test and the estimate 
for the failure rate of the component as in the procedure. 


IF (TY(1).NE.2) THEN 
SUM = 0.0 
DO 50 K = 1, NF(I) 
SUM = SUM + OFT(K) 
50 CONTINUE 
TT(I1) = FLOAT(NC(1)-NF(1I))*OFT(NF(I)) + SUM 
IF (TY(I).EQ.1) THEN 
ELM(I) = FLOAT(NF(1)-1)/TT(I) 
ELSEIF (TY(1).EQ.3) THEN 
ELM(I) = FLOAT(NF(I)-1)/TT(1) 
ENDIF 
ELSEIF (TY(I).EQ.2) THEN 


PRINT 55, M, I 
55 FORMAT(1X,'REP = ',13, ° Comp = ',13,/) 


CALL MLESHAPE(OFT,NC(I),NF(I) ,TOL,1.D0, EBETA(T) ) 
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EBETA(I) = BN(NC(I))*EBETA(T) 


SUM = 0.0 
DO 60 K = 1, NF(T) 

SUM = SUM + OFT(K)**EBETA(T) 
CONTINUE 


TT(I) = FLOAT(NC(I) -NF(I))*OFT(NF(1I))**EBETA(I) + SUM 
ELM(I) = FLOAT(NF(I))/TT(1I) 
ENDIF 


CONTINUE 


Determine the total number of failed test items. 


80 


ISUM = 0 

DO 80 I = 1, NCOMP 
ISUM = ISUM + NF(I) 

CONTINUE 

NFC(M) = ISUM 


Determine the maximum failure rate estimate 
and identify that component. 


90 


ELMAX(M) = 0.0 
KEY = 0 
DO 90 I = 1, NCOMP 
IF (ELM(I).GT.ELMAX(M)) THEN 
ELMAX(M) = ELM(I) 
KEY = I 
ENDIF 
CONTINUE 


Compute the ratios of the failure rate estimates to their maximum. 


100 


DO 100 I = 1, NCOMP 
ER(I) = ELM(I) / ELMAX(M) 
CONTINUE 


Determination of LMU(M) 


110 


SUM = 0.0 
DO 110 I = 1, NCOMP 

SUM = SUM + (ER(I)*TTCI)) 
CONTINUE 


LMU(M) = CHISQD(1-ALPHA, 2*(NFC(M)-NCOMP) ) /(2*SUM) 


Compute estimate of overall reliability RSL(M) for the system. 
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RSE(M) =1.0 
DO 130 J = 1, NCS 
PROD = 1.0 
DOI OeT w= 1, COMP(J,1) 


K = COMP(J,I+1) 
IF (TY(K).EQ.1) THEN 
PROD = PROD*(1 - SURV(TY(K) , LMU(M)*ER(K) , EBETA(K) , UT(K))) 
ELSEIF (TY(K).EQ.2) THEN 
PROD = PROD*(1 - SURV(TY(K) , (LMU(M)*ER(K) )**(1./EBETA(K)) , 
* EBETA(K) , UT(K))) 
ELSEIF (TY(K).EQ.3) THEN 
PROD = PROD*(1 - SURV(TY(K),1.DO-LMU(M)*ER(K) ,0.DO,UT(K))) 
ENDIF 
120 CONTINUE 
REL2(J) = 1.0 - PROD 
RSL(M) = RSL(M) * REL2(J) 
130 CONTINUE 


Increment replication counter. 


SUBROUTINE EVAL 


This subroutine calls BUBBLE to sort the array RSL(NREP) in 
ascending order to get an ordered array ORSL(NREP). It also 
determine the estimate for RSLOW at the specified significance 
level ALPHA and the value of LEVEL in which ORSL(LEVEL) is closest 
to the true reliability RS. 


Include the declaration files 
and declare the local variables. 


INCLUDE 'PARM1 DEF’ 
INTEGER INDEX 
REAL*8 DIFF 
Order the array RSL(NREP) in ascending order. 
DO 10 M = 1, NREP 
ORSL(M) = RSL(M) 
10 CONTINUE 


Bubble Sort. Sink the larger of the pair. 
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CALL BUBBLE(NREP,RSL,ORSL) 


Determine the (1-ALPHA)% lower confidence bound for the system 
reliability. 


RSLOW = ORSL(NINT(NREP*(1-ALPHA) ) ) 


Finding the % confidence level for the true reliability RS. 
(ie. the proportion of RSL(M) lesser than RS) 


DIFF —1,0 
INDEX = 0 
DO 200 M = 1, NREP 
IF (ABS(ORSL(M)-RS).LT.DIFF) THEN 
DIFF = ABS(ORSL(M) -RS) 
INDEX = M 
ENDIF 
200 CONTINUE 


LEVEL = FLOAT( INDEX) /NREP 
Record evaluated parameters in RAW1.DAT (unit 2). 


OPEN (UNIT=2 ) 

WRITE(2, 300) 
300 FORMAT(1X, ' M LMU (M) ELMAX (M) RSL(M) ', 

* : ORSL(M) NFC(M)’ ) 

DO 500 M = 1, NREP 

WRITE(2,400) M,LMU(M) , ELMAX(M) ,RSL(M) ,ORSL(M) , NFC(M) 

400 FORMAT(1X,16,2F12.7,2F12.7,110) 
500 CONTINUE 

CLOSE (UNIT=2 ) 


RETURN 
END 
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SUBROUTINE REPORT 


This subroutine record the simulation results into the ‘'OUT1.DAT’ 
file as logic unit 3. 


Include the declaration files 
and declare local variables. 


INCLUDE 'PARM1 DEF’ 
INTEGER I, J, K, DUM(10) 
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Geewrite to output file ‘OUT1.DAT’ designated as logic unit 3. 
C 

OPEN(UNIT=3 ) 
C 

WRITE(3,10) 

WRITE(3,20) NREP 

WRITE(3,25) 

WRITE(3, 26) 

WRITE(3, 30) 

WRITE(3 ,40) 

WRITE(3,50) ISEED,NCOMP,ALPHA, TOL,NCS , TCN 


WRITE(3,60) 
DO 200 I = 1, NCOMP 
WRITE(3,70) I,TY(1),PARM(1,1) ,PARM(2,1) ,UT(1) ,NC(I) ,NF(I) 
200 CONTINUE 
WRITE(3,80) 
WRITE(3,90) 
DO 300 I = 1, NCOMP 
WRITE(3,100) I,NF(1),TT(1),ELM(I) ,ER(1) , EBETA(I) 
300 CONTINUE 
WRITE(3,110) 
WRITE(3,120) 
DO 500 J = 1, NCS 
DO 400 K = 1, 10 
DUM(K) = COMP(J,K) 
400 CONTINUE 
WRITE(3,130) J,DUM,REL1(J) ,REL2(J) 
500 CONTINUE 
WRITE(3,140) 
WRITE(3,150) RS,ELMAX(NREP) , LMU(NREP) , RSLOW, LEVEL 


C 
10 FORMAT(1X,’OUT1.DAT : Output File of the RETP1 simulation’ ) 
20 FORMAT(1X,’ after ',15,’ replications’ ,/) 
25 FORMAT(1X,’COMMENTS : 8 COMPONENTS IN SERIES ‘) 
26 FORMAT(1X,’ DF = 2 * (NFC - NCOMP) ‘',/) 
30 FORMAT(1X,’Input Parameters:’ ,/) 
40 FORMAT(1X, ’ ISEED NCOMP ALPHA TOL NCS TONE, 7) 


50 FORMAT(IX,F10.1,18,F8.4,F8.5,216,/) 


60 FORMAT(1X,’ I TY(1I) PARM1(1) PARM2(1) UT(I) NC(T) NF(I)',/) 


VOM FORMAT (1X ,13,16,2F9.5,F8.2,218) 
80 FORMAT(1X,/,’Output Parameters for the LAST Replication: ’ ,/) 


90 FORMAT(1X,’ I NF(I) swe) EIM(TI) ERC 


* : EBETA(I)’ ,/) 
HUG  SEORMATCIX,13,16,E16.7,2F14.7,F14.7) 
110 FORMAT(1X,/,'Cut-Set Data:’ ,/) 
120 FORMAT(1X,’ J NUM Component List sas 
* : REL1 Rec an) 
30" FORMAT CIX,13,15,913,2F12.9) 
140 FORMAT(1X,/,’ RS ELMAX (M) LMU(M)’, 


Sil 


* ' RSLOW LEVEL’ , /) 
150 FORMAT(1X,5F12.7,/) 


CLOSE(UNIT=3) 


RETURN 
END 
This portion of the file contains functions and subroutines 
used in the RETP1 model. 
- 14 May 91 
- by Yee Kah-Chee SMC 2802 
A. Random Number Generating Subroutine (LRNDPC). 
(Courtesy of Mr. David Lim Hung-Heng) 
SUBROUTINE LRNDPC (DSEED,U,N) 
INTEGER Niel 
REAL*8 UCN) 
REAL*8 D3SIMT YS DSEED? pot 
D31M1=2**31 - 1 
D31 =2**31 
DATA D31M1/2147483647 .D0/ 
DATA D31 /2147483648.D0/ 
DO 5 I=1,N 
DSEED = DMOD( 950706376 .DO*DSEED,D31M1) 
DSEED = DMOD(16807 .DO*DSEED, D31M1) 
5 U(I) = DSEED / D31 / 
RETURN 
END 


FUNCTION SURV(TYPE, PAR1, PAR2, UTIL) 


This function returns the survival probability of the component of 
different types (TYPE) with scale (PAR1) and shape (PAR2) parameters 
given the specified utilization times or cycles (UTIL). 


INTEGER TYPE, N 
REAL*8 PAR], PAR2, UTIL 


IF (TYPE.EQ.1) THEN 
SURV = EXP(-(PAR1*UTIL) ) 
ELSEIF (TYPE.EQ.2) THEN 
SURV = EXP(-((PAR1*UTIL)**PAR2) ) 
ELSE 
N = NINT(UTIL) 
SURV = PAR1**N 
ENDIF 
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SUBROUTINE BUBBLE(N, LIST, OLIST) 


This subroutine performs a bubble sort in increasing order (ie. sink 
the greater numeral) for the first N terms in an array LIST and 
returns the result in OLIST. 


LOGICAL DONE 
INTEGER N, K, PAIR 
REAL*8 LIST(*), OLIST(*) 


Sink the larger of the pair. 


po 50 K=1, N 
OLIST(K) = LIST(K) 
50 CONTINUE 
PAIR=N-1 
DONE = .FALSE. 
DO WHILE (.NOT.DONE) 
DONE = .TRUE. 
DO 100 K = 1, PAIR 
IF (OLIST(K) .GT.OLIST(K+1)) THEN 
TEMP = OLIST(K) 
Olict (hme OLIST(K+1) 
OLIST(K+1) = TEMP 
DONE = .FALSE. 
ENDIF 
100 CONTINUE 
PAIR = PAIR - 1 
ENDDO 
END 


FUNCTION BN(T) 


This functon returns the value of the unbiased factor for the biased 
maximum likelihood estimate of the shape parameter of a Weibull 
distribution with a sample size of N. 


INTEGER I 
IF (I1.LE.5) THEN 
BN = (1*0.699)/(5.0) 
ELSEIF (1.EQ.6) THEN 
BN = 0.752 
ELSEIF (1.EQ.7) THEN 


so 


BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 


BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSelr 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
EESEEE 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSE 
BN = 
ENDIF 
RETURN 
END 


0.786 

(I.EQ.8) THEN 

0.82 

(I.EQ.9) THEN 

0.8395 

(I.EQ.10) THEN 

0.859 

(I.EQ.11) THEN 

0.871 

(I.EQ.12) THEN 

0.883 

(I.EQ.13) THEN 

0.892 

(I1.EQ.14) THEN 

0.901 

(I1.EQ.15) THEN 

0.9075 

(I1.EQ.16) THEN 

0.914 

(I1.EQ.17) THEN 

0.9185 

(I.EQ.18) THEN 

0.923 

(I.EQ.19) THEN 

0.927 

(I1.EQ.20) THEN 

0.931 

(I1.LE.25) THEN 
0.931+(1-20)*0.014/5.0 
(I.LE.30) THEN 
0.945+(1-25)*0.01/5.0 
(I.LE.40) THEN 
0.955+(1I-30)*0.011/10.0 
(I.LE.60) THEN 
0.966+(I-40)*0.012/20.0 
(I.LE.80) THEN 
0.978+(1-60)*0.006/20.0 
(I.LE.100) THEN 
0.984+(I-80)*0.003/20.0 
(I.LE.120) THEN 
0.987+(I-100)*0.003/20.0 
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SUBROUTINE MLESHAPE(T,N,R,DEL,B, BNEW) 
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This subroutine returns a biased estimator (BNEW) for a Weibull 
shape parameter using the Newton-Raphson’s Method of Successive 
Approximation. The data parameters consist of an ascending ordered 
list of failure times (T), sample size (N), number of failed samples 
(R), tolerance for convergence (DEL) and an initial estimate of the 
shape parameter (B). 


LOGICAL DONE 

INTEGER N, R, I 

REAL*8 GFUNCT, GPRIME, B, BOLD, BNEW, T(*), DEL, 
* TERM1, TERM2, TERM3, SUM1, SUM2, SUM3, SUM4, STEP 


BNEW = B 
DONE = .FALSE. 


DO WHILE (.NOT.DONE) 


DONE = .TRUE. 

TERM] = FLOAT(N-R)*(T(R)**BNEW) 

TERM2 = FLOAT(N-R)*(T(R)**BNEW) *LOG(T(R) ) 

TERM3 = FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R))*LOG(T(R)) 


SUM1 = 0.0 
SUM2 = 0.0 
SUM3 = 0.0 
SUM4 = 0.0 


DO 50 1=1,R 
SUM1 = SUM1 + T(1)**BNEW 
SUM2 = SUM2 + (T(I)**BNEW)*LOG(T(I)) 
SUM3 = SUM3 + (T(I)**BNEW)*LOG(T(1))*LOG(T(T)) 
SUM4 = SUM4 + LOG(T(I)) 
50 CONTINUE 


GFUNCT = (SUM2+TERM2)/(SUM1+TERM1L) - (1.0/BNEW) 


* - (1.0/FLOAT(R) )*SUM4 

GPRIME = (1.0/(SUM1+TERM1)**2)*( (SUM1+TERM1 )*(SUM3+TERM3) 
* - (SUM2+TERM2)**2 ) 
A + (1.0/BNEW**2 ) 


PRINT 60, GFUNCT,GPRIME, BNEW 
60 FORMAT(1X, ‘GFUNCT =’ ,F8.3,'’ GPRIME =’,F8.3,' BNEW =’ ,F8.3) 


Control magnitude of the marching step towards convergence 
as no more than 0.1. 


IF ((GFUNCT.LT.0) .AND. (GPRIME.GT.0)) THEN 


STEP = VMAX(-.1D0, (GFUNCT/GPRIME) ) 
ELSEIF ((GFUNCT.GT.0) .AND. (GPRIME.LT.0)) THEN 
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STEP = VMAX(-.1D0, (GFUNCT/GPRIME) ) 


ELSE 
STEP = VMIN(.1D0, (GFUNCT/GPRIME) ) 
ENDIF 
BOLD = BNEW 
BNEW = BNEW - STEP 


Check for convergence of the MLE for the shape parameter B. 


IF (ABS (BOLD-BNEW) .GT.DEL) THEN 
DONE = .FALSE. 
ENDIF 


Avoid overflow error due to large MLE value caused by small 
GPRIME (slope) as GFUNCT approaches to near zero. 
Stop when magnitude of BNEW exceeds 7. 


IF (BNEW.GT.7.0) THEN 


BNEW = BOLD 
DONE = .TRUE. 
ENDIF 
ENDDO 
RETURN 
END 


eee f®# w®e w®# fF @w®F fF fF fF 2®F fF fe fe 2BF 2B fF @®e weeewprpedqweeeeeteeeee# @e e2eeeweeeeeeeees#e #f8e w@weweeeeeee@qGgeeqeweewweewveeeweeee=2z = 


FUNCTION CHISQD(P,N) 


Modified version of Algorithm 451 from Comunications of the ACM 
Aug 1977 Vol.16 No.8 


This function evaluates the quantile at the probability level P 
(left tail area) for the Chi-square distribution with 
N degrees of freedom. 


REAL*8 P 

REAL X 

INTEGER IF 

DIMENSION C(21), A(19) 

DATA C/ 1.565326E-3, 
1.060438E-3, 
-6.950356E-3, 
“1, 323203h200 
2.277679E-2, 
-8.986007E-3, 
-1.513904E-2, 
2.530010E-3, 


ee * % % 
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10 


20 


30 


a ee 
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+ % s+ % 


el 4S Oluligiee on 
5.169654E-3, 
=f. 153761R <2. 
1.128186E-2, 
2.607083E-2, 
-0.2237368, 

9.780499E-5, 
-8.426812E-4, 
3.125580E-3, 
-8.553069E-3, 
1.348028E-4, 
0.4713941, 

1.0000886 / 


DATA A/ 1.264616E-2, 


-1.425296E-2, 
1.400483E-2, 
-5.886090E-3, 
-1.091214E-2, 
-2.304527E-2, 

Be o4 LIE 3: 
-2.7/28484E-4, 
-9.699681E-3, 

1.316872E-2, 

2.618914E-2, 

BO 22227225 
5.406674E-5, 
3.483789E-5, 

-7.274761E-4, 
3.292181E-3, 

-8.729713E-3, 
0.4714045, 

Ne. 

IF (N-2) 10, 20, 30 

CALL XFROMP(.5*(1.-P) ,X,IF) 

CHISQD = X 

CHISQD CHISQD*CHISQD 

RETURN 

CHISQD 

RETURN 

F = N 

Fl = 1./F 

CALL XFROMP(P,X,IF) 

T =X 

F2 = SQRT(F1)*T 

IF (N.GE. (2+INT(4.*ABS(T)))) GO TO 40 

CHISQD = (((((((C(1)*F2+C(2) )*F24+C(3))*F24+C(4) )*F2 
+C(5))*F2+C(6) )*F24+C(7) )*F1+(( (CC (C(8)4+C(9) *F2) *F2 
+C(10) )*F2+C (11) )*F2+C(12) )*F2+C(13))*F2+C(14) ) )*F1+ 
(CCC (C(15) *F24+C (16) )*F2+C(17) )*F2+C(18) ) *F2 
+C(19) )*F24+C(20) )*F2+C(21) 


=) SOG Cle +P) 


ay 
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GO TO 50 
40 CHISQD=(((A(1)+A(2)*F2)*F1+(((A(3)+A(4) *F2) *F2 
+A(5))*F2+A(6)))*F1+(((((A(7)+A(8)*F2)*F2+A(9) )*F2 
+A(10))*F2+A (11) )*F2+A(12)))*F1+(((((A(13) *F2 
+A(14))*F2+A(15) )*F2+A(16))*F2+A(17) )*F2*F2 
+A(18))*F2+A(19) 
50 CHISQD = CHISQD*CHISQD*CHISQD*F 

RETURN 

END 
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SUBROUTINE XFROMP(P,X,IFAULT) 
Algorithm AS 24 J.R.STAT.SOC. C. (1969) Vol.18. No.3. 


This subroutine computes the standard normal deviate X for 
the specified left tail area P. 


REAL*8 P 

DIMENSION A(5) 

DIMENSION CONNOR (17), HSTNGS(6) 
DATA CONNOR/ 8.0327350124E-17, 
.4483264644E-15, 
.4668270103E-14, 
-9554295164E-13, 
.9477940136E-12, 
.3507027951E-11, 
.0892221037E-9, 
.3122532964E-8, 
.4503852223E-7, 
.4589169001E-6, 
.3227513228E-5, 
.0683760684E-4, 
.5757575758E-4, 
.6296296296E-3, 
.3809523810E-2, 
ig Ls 

33333333333 / 
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DATA RTHFPI / 1.2533141373 / 
DATA RRT2PI / 0.3989422804 / 
DATA TERMIN / 1.0E-11 / 


DATA HSTNGS / 2.515517, 


* 02602853. 
* 0.010328, 
* Pas 2 7 oo 
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* 0.189269, 
* 0.001308 / 


IFAULT = 1 
IF ((P.LE.0.0).OR.(P.GE.1.0)) GO TO 100 
IFAULT = 0 


first approximation XO to deviate by Hastings’ formula 


B = P 
IF(B.GT.0.5) B= 1.0 - B 


F = - LOG(B) 

E = SQRT(F+F) 

XO = -E + ((HSTNGS(3)*E+HSTNGS (2) )*E+HSTNGS(1))/ 
(( (HSTNGS (6) *E+HSTNGS (5) )*E+HSTNGS (4) ) *E+1 .0) 

i (XOLLT.0.0) COnTO. 1 

XO = 0.0 

PO = 0.5 

X1 = -RTHFPI 

iO ae) 7 


Find the area PO corresponding to XO 


1 


(1) 


(2) 


Y = XO%x?2 
IF (XO.LE.-1.9) GO TO 3 
Y = -0.5*Y 


series approximation 


PO = CONNOR(1) 
DO 2 L=2,17 
PO = PO*Y + CONNOR (L) 


PO = (PO*Y+1.0)*XO 

Xl = -(PO+RTHFPI)*EXP(-Y) 
PO = PO*XRRT2PI + 0.5 
GOrro! 7 


continued fraction approximation 


Z=1.0/Y 

A(2) = 1.0 
A(3) = 1.0 
A(4) = Z+1.0 
A() =:1.0 

W = 2.0 

DO 6 L=1,3,2 
DO 5 J-1,2 
hae, 


Sy 


KA = 7 - K 
5 ACK) = A(KA) + A(K)*W*Z 


6 W=W+ 1.0 
APPRXU = A(2)/A(3) 
APPRXL = A(5)/A(4) 
C = APPRXU - APPRXL 
IF (C.GE.TERMIN) GO TO 4 
X1 = APPRXL/XO 
PO = -X1*RRT2PI*EXP(-0.5*Y) 
C 
C Get accurate value of deviate by Taylor Series 
C (X1, X2, X3 are derivatives for the Taylor Series 
C 
7 D = F + LOG(PO) 
X2 = XO*X1*X1 -X1 
X3 = X1¥**3 + 2.0*XO*X1*X2 -X2 
X = ((X3*D/3.0+X2)*D/2.0+X1)*D + XO 
IF (P.LE.0.5) GO TO 100 
X = -X 
100 RETURN 
END 


FUNCTION VMAX(X, Y) 
REAL*8 X, Y 
IF(X.GT.Y) THEN 
VMAX = X 
ELSE 
VMAX = Y 
ENDIF 
RETURN 
END 


FUNCTION VMIN(X,Y) 

REAL*8 X,Y 

LRCx2 iY) THEN 
VMIN = X 

ELSE 
VMIN = Y 

ENDIF 

RETURN 

END 


4. Program Output. (OUT1.DAT) 
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The result for the simulation run based on the input parameters specified in 
IN1.DAT are computed and written to the file OUT1.DAT. A sample of this file is 
as follows. 


OUT1.DAT : Output File of the RETP1 simulation 
after 1000 replications 
COMMENTS : 8 COMPONENTS IN SERIES 


DF = 2 * (NFC - NCOMP) 


Input Parameters: 
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ISEED NCOMP ALPHA TOL NCS TCN 
16807 .0 8 0.2000 0.01000 8 

I TY(I) PARM1(1I) PARM2(T) UT(I) NC(I) NF(IL) 

i 1 0.00200 1.00000 10.00 1) 3 

2 1 0.00400 1.00000 10.00 iis: 3 

3 1 0.00600 1.00000 10.00 i | 

4 1 0.00800 1.00000 10.00 15 3 

5 2 0.00200 2.00000 10.00 15 3 

6 2 0.00400 2.00000 10.00 15 3 

7 2 0.00600 2.00000 10.00 ales) 3 

8 2 0.00800 2.00000 10.00 iS 3 

Output Parameters for the LAST Replication: 
I NF(I) ETC ELM(1) ER(1) EBETA(L1) 
I S) 0.1969091E+04 0.0010157 0.1703872 O.0000000 
2 3 0.3355078E+03 0.0059611 1 .0000000 0 .0000000 
3 3 0.5999981E+03 0 .0033333 OSS ovale ls. 0 .0000000 
4 3 0.6423594E+03 0.0031135 Oe 5223055 0 .0000000 
5 3 0.2563836E+16 0 .0000000 0 .0000000 6.3524983 
6 ) 0.8186629E+07 0 .0000004 0.0000615 22957) loaG 
7 3 02367 5037E+09 0 .0000000 0 .0000009 4.2418658 
8 3 Ooo sD 19EFOS 0.0000529 0.0088672 2.3471863 
Cut-Set Data: 

J NUM Component List RELL REL2 (M) 

1 li 00 0 0 0 0 0 O 0.980198622 0.990279973 

2 120 0 0 0 0 0 0 O 0.960789382 0.944286704 

3 13 000 0 0 0 0 O 0.941764534 0.968452990 

4 On) 0 OOO mmOmmOnO 9231632600 ,970502377 

5 1500 0 0 0 0 0 O 0.999600053 0.999999940 

6 16000 0 0 0 0 O 0.998401225 0.999680758 


170 00 0 0 0 0 O 0.996406436 0.999911249 
18 0 0 0 © 0 0 0 © 0.99362037%6rt- 3287373 
RS ELMAX (M) LMU (M) RS LOW LEVEL 


0.8089645 0.0059611 0.0057325 0.8691733 0.3380000 
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APPENDIX C: Users’ Guide for RETP2 


Reliability Estimation Test Plan 2 (RETP2). 
by YEE, Kah-Chee 


July 91 





1. Brief Description. 


RETP2 is a computer program written in FORTRAN that runs on the Amdahl 
mainframe at NPGS. It allows the user to simulate exponential and Weibull failure 
times of component items being tested to evaluate the accuracy of a confidence limit 
estimation procedure based on Type I data censoring (that is, testing items of 
component i until a specified total testing time is achieved for each of them). 


2. Program Input. (IN2.DATA) 


The input of the program are specified to the program via an input file called 
IN2.DAT. A sample input file is shown below. 


This file contains the inputs required by the RETP2 model. 
Update only the numerical values between dotted lines as appropriate. 
Do not delete any of the comment lines. (IN2.DAT) 20 Jun 91 


C Value Type Units Description Variable 
oso o notes. B boise 5 oes OS es 
16807 .0 Real - initial random seed ISEED 
8 Int - total # of components in system NCOMP 
4 INT - j## OF EXPONENTIAL COMPONENTS NEXP 
4 INT - ## OF WEIBULL COMPONENTS NWEI 
0 INT - ## OF GEOMETRIC COMPONENTS NGEO 
0.01 Real - tolerance for MLE TOL 
0.20 REAL - DESIRED SIGNIFICANCE LEVEL ALPHA 
1000 INT - ## OF REPLICATIONS DESIRED NREP 
5 INT - TEST CASE NUMBER TCN 
1 = all EXP 
2 = all WEI 


3 = EXP + WEI 
4 = EXP + WEI + CYC 
8 Int - number of cut sets NCS 


C TEST PLAN : Testing until TT(1) (total test time) is accumulated. 
C (Use REAL numbers ONLY!!!) 
Gow ww wwe tees en eee eb eee ce eee ee ee ea 6 6S eee see ee ee es ee ee 
C Comp Comp Comp Parameters Util Test Plan Inputs 
C Number Type Scale Shape Time/Cycle Accumulated 
Cc 61 TY(I) PARM1(I1) PARM2 (1) UT(I) i NC(I) 
C Ime nt Real Real (hrs) (cycs) (hrs) (cycs) Int 
C--4-- 0 oo oe ee eee ee tee See = eee. 2 eee 
130 ee) 0.005 lee) Db. 0 5400.0 20.0 
2.0 be) 0.005 1.0 5.0 5400.0 20.0 
3.0 1.0 0.005 lus’ 5.0 5400 .0 20.0 
4.0 iy..0 0.005 0 5.0 5400.0 20.0 
5.0 2.0 0.010 2.0 L5igO 2700.0 20.0 
6.0 2.0 0.010 2.0 15.0 2700.0 20.0 
7.0 2.0 0.010 2.0 1320 2700.0 20.0 
8.0 2.0 0.010 200) 15.0 2700.0 20.0 
Cn on nn i ee oe 
C Note : TY(I)=l EXPONENTIAL P(surv) = exp(-PARM2)*T) 
G TY(I)=2 WEIBULL P(surv) = exp(- (PARM1*T)**PARM2) 
C TY(I)=3 GEOMETRIC P(surv) = PARM1**T 
oe ee ese ee ec 
C SYSTEM CONFIGURATION : Identification of CUT SETs 
C - min groups of components that have to fail 
C for the system to fail. 
nn ene ee ee 
C Cut Set # in Set List of Components in Cutset 
C J COMP(J,1) COMP(J,2) ... up to COMP(J,1) components 
Cann ne ee eee we ee ne ee es oo see eee 
i 1 1 0° 0° "OO 0 OTH Or @ 
2 1 2000 0 0 0 0 0 0 
3 J. 3 0 00 0 0 0 0 0 0 
4 i 4 00 0 0 0 0 0 0 0 
5 il 5 0 00 00 0 0 0 0 
6 an 6 0 0 0. 0-0 FO 10m OD 
7 1 7 0 0 0 0 0 0 0 0 0 
8 Jk 8 0 0 C0 0 0 0 0) OMG 
Conce eee e ee eee ce cn ee ee eee 


3. Program Flow and Logic. (NAME2.DEF, PARM2.DEF and RETP2.FOR) 


Input parameters are first read in by the program by calling the INPUT 
subroutine. The program then evoke the SIM subroutine which generates the 
random failure times and compute the key statistics required in the procedure. The 
next subroutine EVAL determines the measures of accuracy for run. REPORT is 
the subroutine which generates the output file for the run OUT2.DAT. 


The variables in the program RETP2.FOR are described in the file 
NAME2.DEF as listed below. 
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This file contains the declaration for input and output variables 
used in the the RETP2 model. (NAME2.DEF) 20 Jun 91 


77s eweeegweeeeFweeeweeweeeeweerweweewe ww @ eevee @f®@ we w#ewFeweeweeweweeweweewewe w# ww we we eweeeweeweeweeweeweweweeweweweeweweeeweewnsweeeweweseeese see @ = = 


ISEED = initial random seed selected. 

SEED = current random seed. 

RS = true overall series system reliability. 

ALPHA = level of significance desired. 

NREP = number of replications desired for the simulation. 
TPN = test plan number (1). 

TCN = test case number (l, 2, 3 or 4). 

NCOMP = total number of components in the system. 

NEXP = number of components with EXP failure times. 
NWEI = number of components with WEI failure times. 
NGEO = number of components with GEO failure times. 
TOL = desired tolerance for MLE of WEI shape parameter. 


Distribution: EXPonential WEIbull GEOmetric 

TY(I) = type: 1 2 3 
PARM(1,1): Scale(l/hr) Scale(l/hr) Prob 
PARM(2,I1): - Shape - 


UT(I) = utilization time (hrs) for component i (EXP and WEI). 

UC(I) = utilization cycles for component i (GEO only). 

NC(I) = number of test samples (sample size) for cmponent i. 

NF(I) = desired number of failures in test for component i. 

NCS = number of cut-sets for the system. 

COMP(J,K) = kth parameter of cut-set j (first being the no. of 
components belonging to the cut-set). 


Assumed Variables. 

MAXCOMP = maximum number of components allowed in the system. 
MAXREP = maximum number of replications permitted. 

MAXCUT = maximum number of cut-sets. 


Program and Output Variables. 
RS = true overall system reliability. 
TT(I) = total accumulated failure time (hr) for component i 
(EXP and WEI only). 
TC(I) = total accumulated cycles to failure (incl. failure cycle) 
for component i (GEO only). 
EBETA(I) = estimate for shape parameter of component i (if Weibull). 
REL1(J) = actual reliability for cut-set j. 
REL2(J) = computed reliability for cut-set j for current replication. 
ELM(I) = estimated component failure rate (l/hrs) for component i. 
ELMAX(M) = max estimated component failure rate for rep m (l/hrs). 
ER(I) = ratio of estimated failure rate to LMAX. 
ET(I1) = same as TT(I1) except that these are for Weibull components. 
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NFC(M) total number of failed test components. 

LMU (M) upper confidence bound for failure rate (l/hrs), 

RSL(M) = lower confidence limit estimated for system reliability 

for the mth replication. 

ORSL(M) = ordered RSL(M) (ascending). 

RSLOW = (1-ALPHA)x100 percentile of set of RSL(M). 

LEVEL = achieved confidence level, ie. proportion of RSL(M) that 
are lesser than RS (conservative estimate). 


-- END OF NAME2.DEF == -- ~~ ~~ 9 ge( ge cps oP icgs ieee = ae 
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Together with the main program in RETP2.FOR are the other subroutines 
needed in the simulation. The declaration of variables is done in the file 
PARM2.DEF. Relevant descriptions are included as comment lines in the source 


code to help explain the program segments. A listing of PARM2.DEF and 
RETP2.FOR is given below. 


C This file contains the declaration for input and output variables 
C used in the the RETP2 model. (PARM2.DEF) 20 Jun 91 


INTEGER MAXCOMP, MAXREP 

PARAMETER( MAXCOMP = 100 , MAXREP = 1000 , MAXCUT = 20 ) 

REAL*8 ISEED, SEED 

INTEGER NREP, TCN, NCOMP, NEXP, NWEI, NGEO, NCS, 
NC(MAXCOMP) , NF(MAXCOMP), TY(MAXCOMP), NFC(MAXREP) , 
UC(MAXCOMP) , TC(MAXCOMP), COMP(MAXCUT ,MAXCOMP) 

REAL*8 RS, ALPHA, UT(MAXCOMP), TT(MAXCOMP) , 


+ % 


* PARM(2 ,MAXCOMP), ELM(MAXCOMP), ER(MAXCOMP), ET (MAXCOMP) , 

* LMU (MAXREP), RSL(MAXREP), ORSL(MAXREP) , 

* ELMAX (MAXREP) , RSLOW, LEVEL, TOL, EBETA(MAXCOMP) , 

* REL1(MAXCUT) , REL2 (MAXCUT) 
C 

COMMON/BLOCK1/ISEED, SEED, NREP, TCN, NCOMP, NC, NF, NEXP, NWEI, 
* NGEO, NCS, TY, NFC, UC, TC, COMP 
COMMON/BLOCK2/RS, ALPHA, UT, TT, PARM, ELM, ER, ET, LMU, 

* RSL, ORSL, ELMAX, RSLOW, LEVEL, TOL, EBETA, 

* REL1, REL2 
C 
C-- END OF PARM2.DEF ---------------------- rrr rr rr rrr rr rr rrr crn 
C 
(ee ee 
C This file contains the main program and the subroutines 
C for the Reliability Estimation Test Plan 2 (RETP2) model. 
C (RETP2.FOR) - runs on a IBM PC Compatible. 
C 
C IBM Mainframe Version. 
C 
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Test Plan 2 : Testing until accumulated time or cycles is achieved 


----------- for component i. 


by Yee Kah-Chee SMC 2802. 
20 Jun 91. 


PROGRAM RETP2 
Include the declaration files. 


INCLUDE ‘'NAME2 DEF’ 
INCLUDE 'PARM2 DEF’ 


Read in input data. 
CALL INPUT 
Activate simulation. 
CALL SIM 
Process and evaluate output data. 
CALL EVAL 
Generate simulation report. 
CALL REPORT 


STOP 
END 
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SUBROUTINE INPUT 


This subroutine reads in the inputs for the RETP2 model. 


Include the declaration file. 
INCLUDE ‘'PARM2 DEF’ 


INTEGER I, J, K, DUM2(11) 
REAL*8 DUM1(7) 


Read data from 'IN2.DAT’ designated as logic unit l. 


OPEN(UNIT=1) 
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READ(1,10) 

10 FORMAT(1X,/////) 
READ(1,*) ISEED 
READ(1,*) NCOMP 
READ(1,*) NEXP 
READ(1,*) NWEI 
READ(1,*) NGEO 
READ(1,*) TOL 
READ(1,*) ALPHA 
READ(1,*) NREP 
READ(1,*) TCN 
READ(1, 20) 

20 FORMAT(1X,///) 
READ(1,*) NCS 


READ (1,30) 
30 FORMAT(1X,////////) 


DO 50 K = 1, NCOMP 
READ(1,*) DUM1 


I = NINT(DUM1(1)) 
TY(I) = NINT(DUM1(2) ) 


IF (TY(I).EQ.1) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 


UT(I) = DUM1(5) 
TT(I) = DUM1(6) 
NC(I) = NINT(DUM1(7)) 


EBETA(I) = 0 

ELSEIF (TY(1).EQ.2) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 
UT(I) = DUM1(5) 
TT(I) = DUM1(6) 

NC(I) = NINT(DUM1(7)) 
ELSEIF (TY(I).EQ.3) THEN 
PARM(1,1) = DUM1(3) 
PARM(2,I) = DUM1(4) 
UC(I) = NINT(DUM1(5)) 

UT(I) = DUM1(5) 


TT(I) = DUM1(6) 
NC(I) = NINT(DUM1(7)) 
ENDIF 


50 CONTINUE 


READ (1, 60) 
60 FORMAT(1X,///////////) 
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DO 80 I = 1, NCS 
READ(1,*) DUM2 
J = DUM2(1) 
COMP(J,1) = DUM2(2) 
DO 70 K = 1, COMP(J,1) 
COMP(J,K+1) = DUM2(K+2) 
70 CONTINUE 
80 CONTINUE 
CLOSE(UNIT=1) 
RETURN 
END 


SUBROUTINE SIM 


This subroutine simulates NREP possible outcomes of the test plan 
desired in order to obtain the raw estimates of LMU(M) and RSL(M) 
for each of the replication. 


Include the declaration file 
and declare local variables. 


INCLUDE 'PARM2 DEF’ 
INTEGER I, J, K, M, ISUM, KEY, ICOUNT 
INTEGER NCYC(MAXCOMP), NSYS 


REAL*8 UNI 

REAL*8 SUM, PROD, 

* FT (MAXCOMP,200), OFT(MAXCOMP, 200), 
* DFT(200), DOFT(200) 


LOGICAL FLAG 
SEED = ISEED 
Compute overall true system reliability RS. 


RS = 1.0 
DO 30 J = 1, NCS 
PROD = 1.0 
DO 20 I = 1, COMP(J,1) 
K = COMP(J,I+1) 
PROD = PROD*( 1 - SURV(TY(K) ,PARM(1,K) ,PARM(2,K),UT(K)) ) 
20 CONTINUE 
REL1(J) = 1.0 - PROD 
RS = RS * REL1(J) 
30 CONTINUE 


Start of Simulation. 
(Initialize replication counter M). 
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M=1 
DO WHILE (M.LE.NREP) 


C 
C Test Plan : Sample and determine unknown NF(I) 
C --------- (with known NC(I), Weibull case) until TT(I) is reached. 
C 
C Generate NC(I) failure times, put them in ascending order 
C with the smallest failure time on the top of the list. 
C 
DO 60 I = 1, NCOMP 
C 
C Exponential components. Test each component until it fails, check 
C to see if TT(I) is exceeded, if not, carry on testing. 
C 
IF(TY(I).EQ.1) THEN 
SUM = 0.0 
L=1 
DO WHILE (SUM.LE.TT(I)) 
CALL LRNDPC(SEED,UNI,1) 
FT(I,L) = -LOG(UNI) /PARM(1,I) 
SUM = SUM + FT(I,L) 
L=L+t+l1 
ENDDO 
NF(I) = L - 2 
C 
C Weibull components. Generate NC(I) failure times, put them in 
C an ascending order with the smallest failure time on top and 
C determine NF(I). 
C 


ELSEIF(TY(1I).EQ.2) THEN 
DO 40 K = 1, NC(I) 
CALL LRNDPC(SEED,UNI,1) 
FT(I,K) = (1.0/PARM(1,1I))*(-LOG(UNI) )**(1.0/PARM(2,1)) 
DFT(K) = FT(I,K) 
40 CONTINUE 
CALL BUBBLE(NC(I) , DFT,DOFT) 
DO 42 K = 1, NC(I) 
OFT(I,K) = DOFT(K) 
42 CONTINUE 
SUM = 0.0 
NF(I) = 0 
DO 45 K = 1, NC(I) 
IF (OFT(I,K).LT.TT(1I)) THEN 
NF(I) = NF(I) + 1 


ENDIF 
45 CONTINUE 
C 
C GEOMETRIC COMPONENTS. DETERMINE NF(TI). 
C 


ELSEIF(TY(1I) .EQ.3) THEN 
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NF(I) = 0 
DO 50 K = 1, TT(I) 
CALL LRNDPC(SEED,UNI,1) 
IF (UNI.GT.PARM(1,1)) THEN 
NF(I) = NF(I) + 1 
ENDIF 
50 CONTINUE 


ENDIF 
60 CONTINUE 


Determine total number of failed test components as well as 
checking for zero component failure or just failures from a 
SINGLE component (FLAG will be set to .TRUE. if so). 


ISUM = 0 
ICOUNT = 0 
DO 70 I = 1, NCOMP 
IF (NF(1).GT.0) THEN 
ICOUNT = ICOUNT + 1 
ENDIF 
ISUM = ISUM + NF(I) 
70 CONTINUE 


NFC(M) = ISUM 
IF (ICOUNT.LE.1) THEN 
FLAG = .TRUE. 
ELSE 
FLAG = .FALSE. 
ENDIF 


Case A : More than ONE component type experienced failures 
------ in the test. 


IF (.NOT.FLAG) THEN 
Estimate the failure rate of each component. 
DO 90 I = 1, NCOMP 
Exponential components. 
IF (TY(1).EQ.1) THEN 
ET(I) = TT(I) 
ELM(I) = FLOAT(NF(1))/TT(1) 


Weibull components. 


ELSEIF (TY(1).EQ.2) THEN 
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DO 75 K = 1, NC(I) 
DOFT(K) = OFT(I,K) 
75 CONTINUE 
CALL MLESHAPE(DOFT,NC(I) , IMAX(1,NF(1)),TOL,1.D0,EBETA(1I) ) 
EBETA(I) = BN(NC(I))*EBETA(L) 
SUM = 0.0 
DO 80 K = 1, NF(I) 
SUM = SUM + OFT(I,K)**EBETA(1) 
80 CONTINUE 
ET(1) = FLOAT(NC(I)-NF(1I))*TT(1)**EBETA(1) + SUM 
ELM(1) = FLOAT(NF(1))/ET(I) 


Geometric components. 


ELSEIF (TY(1I).EQ.3) THEN 
ET(1) = TT(1) 
ELM(I) = FLOAT(NF(I))/TT(I) 


ENDIF 
90 CONTINUE 


Determine the maximum failure rate estimate 
and identify that component. 


ELMAX(M) = 0.0 
KEY = 0 
DO 100 I = 1, NCOMP 
IF (ELM(1I).GT.ELMAX(M)) THEN 
ELMAX(M) = ELM(I) 
KEY = I 
ENDIF 
100 CONTINUE 


Compute the ratios of the failure rate estimates to their maximun. 
DO 110 I = 1, NCOMP 
ER(I) = ELM(I) / ELMAX(M) 
110 CONTINUE 
Determination of LMU(M) 
SUM = 0.0 
DO 120 I = 1, NCOMP 
SUM = SUM + (ER(I1)*ET(1)) 
P20 CONTINUE 
LMU(M) = CHISQD(1-ALPHA,NINT(1.3*2*(1+NFC(M) )))/(2*SUM) 


Compute estimate of overall reliability RSL(M) for the system. 
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RSL(M) = 1.0 
DO 140 J = 1, NCS 
PROD = 1.0 


DO 130 I = 1, COMP(J,1) 
K = COMP(J,I+1) 
IF (TY(K).EQ.1) THEN 
PROD = PROD*(1 - SURV(TY(K) , LMU(M)*ER(K) , EBETA(K) , UT(K) )) 
ELSEIF (TY(K).EQ.2) THEN 
PROD = PROD*(1 - SURV(TY(K) , (LMU(M)*ER(K) )**(1./EBETA(K)), 
* EBETA(K) ,UT(K))) 
ELSEIF (TY(K).EQ.3) THEN 
PROD = PROD*(1 - SURV(TY(K) ,1.DO-LMU(M)*ER(K) ,0.DO,UT(K))) 
ENDIF 
130 CONTINUE 
REL2(J) = 1.0 - PROD 
RSL(M) = RSL(M) * REL2(J) 
140 CONTINUE 


C 
ENDIF 
C 
C Case B : Where there are at most ONE component which experienced 
C ------ failures during the test. 
C 
IF (FLAG) THEN 
C 
C Determine number of complete systems (NSYS) implied by the test 
C based on cut-set information after first determining the number 
C of mission cycles tested for each component (to the nearest 
C integer) (NCYC(I)). 
C 


ISUM = 0 
DO 150 I = 1, NCOMP 
Nene Clje—s INDCTECL)/ULCL)) 
ISUM = ISUM + NCYC(I) 
ELM(I) = 0 
ER(I) = 0 
EBETA(I) = 1.0 
150 CONTINUE 
ELMAX(M) = 0 
LMU(M) = 0 


NSYS = ISUM 
ISUM = 0 
Dom 70a 
ISUM = 
DO 160 
ISUM 
160 CONTINUE 
ISUM = ISUM / COMP(J,1) 


1, NCS 


= 1, COMP(J,1) 
ISUM + NCYC(COMP(J ,I+1)) 


| AO I 
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NSYS = IMIN(NSYS,ISUM) 
170 CONTINUE 
CG 
For zero failures in the test. 


Qe 


IF (ICOUNT.EQ.0) THEN 
RSL(M) = ALPHA**(1.0/FLOAT(NSYS) ) 


Q 


For failures experienced by a particular component type. 
ELSE 
CALL GETP(NSYS , ALPHA, RSL(M) ) 
ENDIF 
ENDIF 


C Increment replication counter. 


O 


SUBROUTINE EVAL 


This subroutine calls BUBBLE to sort the array RSL(NREP) in 
ascending order to get an ordered array ORSL(NREP). It also 
determine the estimate for RSLOW at the specified significance 
level ALPHA and the value of LEVEL in which ORSL(LEVEL) is closest 
to the true reliability RS. 


Include the declaration files 
and declare the local variables. 


Ory 2G) 0) GG) Or © 


INCLUDE 'PARM2 DEF’ 
INTEGER INDEX 
REAL*8 DIFF 


©) 


C Order the array RSL(NREP) in ascending order. 


DO 10 M = 1, NREP 
ORSL(M) = RSL(M) 
10 CONTINUE 


©) 


Bubble Sort. Sink the larger of the pair. 


CALL BUBBLE(NREP,RSL,ORSL) 
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Determine the (1-ALPHA)% lower confidence bound for the system 
reliability. 


RSLOW = ORSL(NINT(NREP*(1-ALPHA) ) ) 


Finding the % confidence level for the true reliability RS. 
(ie. the proportion of RSL(M) lesser than RS) 


DIFF = 1.0 
INDEX = 0 
DO 200 M = 1, NREP 
IF (ABS(ORSL(M)-RS).LT.DIFF) THEN 
DIFF = ABS(ORSL(M) -RS) 
INDEX = M 
ENDIF 


200 CONTINUE 


LEVEL = FLOAT( INDEX) /NREP 
Record evaluated parameters in RAW2.DAT (unit 2). 


OPEN (UNIT=2) 
WRITE(2, 300) 


300 FORMAT(1X, ' M LMU (M) ELMAX (M) RoLCM)’ , 


* ' ORSL(M) NEGCM)) 
DO 500 M = 1, NREP 
WRITE(2,400) M,LMU(M) , ELMAX(M) ,RSL(M) ,ORSL(M) ,NFC(M) 


400 BORMAT Gites 6, 2612, 7,Z2Fl2./7,110) 
500 CONTINUE 


CLOSE(UNIT=2) 


RETURN 
END 


SUBROUTINE REPORT 


This subroutine record the simulation results into the ‘'OUT2.DAT' 
file as logic unit 3. 


Include the declaration files 
and declare local variables. 


INCLUDE ‘PARM2 DEF’ 
INTEGER I, J, K, DUM(10) 


Write to output file ‘OUT2.DAT’ designated as logic unit 3. 


OPEN (UNIT=3 ) 


75 


200 


300 


400 


500 


100 
110 
120 


130 
140 


150 


WRITE(3,10) 

WRITE(3,20) NREP 

WRITE(3, 25) 

WRITE(3 426) 

WRITE(3, 30) 

WRITE(3,40) 

WRITE(3,50) ISEED,NCOMP ,ALPHA, TOL,NCS , TCN 


WRITE(3, 60) 
DO 200 I = 1, NCOMP 
WRITE(3,70) I,TY(I),PARM(1,I), PARM(2,1),UT(I) , TT(1) ,NC(I) ,NF(I) 
CONTINUE 
WRITE(3, 80) 
WRITE(3,90) 
DO 300 I = 1, NCOMP 
WRITE(3,100) I,NF(1I),ET(I) ,ELM(I) , ER(1) , EBETA(I) 
CONTINUE 
WRITE(3,110) 
WRITE(3,120) 
DO 500 J = 1, NCS 
DO 400 K = 1, 10 
DUM(K) = COMP(J,K) 
CONTINUE 
WRITE(3,130) J,DUM,REL1(J) ,REL2(J) 
CONTINUE 
WRITE(3,140) 
WRITE(3,150) RS,ELMAX(NREP) , LMU(NREP) , RSLOW, LEVEL 


FORMAT (1X,’OUT2.DAT : Output File of the RETP2 simulation’ ) 


FORMAT (1X, ' after ',15,' replications’ ,/) 
FORMAT (1X, ‘COMMENTS : 8 COMPONENT IN SERIES ee 
FORMAT (1X, ' DF = NINT (1.3 * 2 * (1 + NFC)) ‘PD 
FORMAT(1X,'Input Parameters: ’ ,/) 

FORMAT (1X, ’ ISEED NCOMP ALPHA TOL NCS TCN’ ,/) 
FORMAT(1X,F10.1,18,F8.4,F8.5,216,/) 


FORMAT(1X,' I TY(I) PARMI1(I) PARM2(1) UT(I) TG) NCC 


. oe NEC ey) 


FORMAT (1X, 135,16, 2F9%5, Zhe. 27 20a) 
FORMAT(1X,/,'Output Parameters for the LAST Replication: '’ ,/) 
FORMAT(1X,’' I NFCT) EECL) ELM(T) ER(CI)’, 


* EBETA(I)’,/) 


* 


* 


FORMAT (1X,1I3,16,E16.7,2F14.7,F14.7) 
FORMAT(1X,/,'Cut-Set Data:’,/) 
FORMAT(1X,' J NUM Component List - 


' REL1 REL2(M)',/) 
FORMAT (1X,13,15,913,2F12.9) 
FORMAT (1X, /, ’ RS ELMAX (M) LMU(M)', 
' RSLOW LEVEL’ , /) 


PORMAT CR Or UZ: fa) 
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CLOSE(UNIT=3 ) 


RETURN 
END 


This portion of the file contains functions and subroutines 
used in the RETP2 model. 
- 20 Jun 91 
- by Yee Kah-Chee SMC 2802 
A. Random Number Generating Subroutine (LRNDPC). 
(Courtesy of Mr. David Lim Hung-Heng) 
SUBROUTINE LRNDPC (DSEED,U,N) 
INTEGER N, I 
REAL*8 UCN) 
REAL*8 D31M1, DSEED, D31 
D31M1=2**31 - 1 
D319 =2**31 


DATA D31M1/2147483647 .D0/ 

DATA D31 /2147483648.D0/ 

DO 5 I=1,N 
DSEED = DMOD(950706376.DO*DSEED, D31M1) 
DSEED = DMOD(16807.DO*DSEED,D31M1) 

5 U(1) = DSEED / D31 
RETURN 
END 
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FUNCTION SURV(TYPE, PAR1, PAR2, UTIL) 


This function returns the survival probability of the component of 
different types (TYPE) with scale (PAR1) and shape (PAR2) parameters 
given the specified utilization times or cycles (UTIL). 


INTEGER TYPE, N 
REAL*8 PAR1, PAR2, UTIL 


IF (TYPE.EQ.1) THEN 
SURV = EXP(-(PARI1*UTIL) ) 
ELSEIF (TYPE.EQ.2) THEN 
SURV = EXP(-((PAR1*UTIL)**PAR2) ) 
ELSE 
N = NINT(UTIL) 
SURV = PAR1**N 
ENDIF 


END 
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SUBROUTINE BUBBLE(N,LIST,OLIST) 


This subroutine performs a bubble sort in increasing order (ie. sink 
the greater numeral) for the first N terms in an array LIST and 
returns the result in OLIST. 


LOGICAL DONE 
INTEGER N, K, PAIR 
REAL*8 LIST(*), OLIST(*) 


Sink the larger of the pair. 


DO 50 K= 1, N 
OLIST(K) = LIST(K) 
50 CONTINUE 
PAIR=N- 1 
DONE = .FALSE. 
DO WHILE (.NOT.DONE) 
DONE = .TRUE. 
DO 100 K = 1, PAIR 
IF (OLIST(K).GT.OLIST(K+1)) THEN 
TEMP = OLIST(K) 
OLIST(K) = OLIST(K+1) 
OLIST(K+1) = TEMP 
DONE = .FALSE. 
ENDIF 
100 CONTINUE 
PAIR = PAIR - 1 
ENDDO 
END 


FUNCTION BN(1) 


This functon returns the value of the unbiased factor for the biased 
maximum likelihood estimate of the shape parameter of a Weibull 
distribution with a sample size of N. 


INTEGER I 
IF (1.LE.5) THEN 
BN = (I1*0.699)/(5.0) 
ELSEIF (1.EQ.6) THEN 
BN = 0.752 
ELSEIF (1.EQ.7) THEN 
BN = 0.786 
ELSEIF (1.EQ.8) THEN 
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BN = 
BLESEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSEIF 
BN = 
ELSE 
BN = 
ENDIF 
RETURN 
END 
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(I.EQ.9) THEN 


0.8395 
(I.EQ. 
0.859 
(1.EQ. 
0.871 
(I.EQ. 
0.883 
(I.EQ. 
0.892 
(1.EQ. 
0.901 
(Je). 
0.9075 
(I.EQ. 
0.914 
(I.EQ. 
0.9185 
(1.EQ. 
0.923 
(I.EQ. 
0.927 
(1.EQ. 
0.931 
‘Cie 


0.931+(1-20)*0.014/5.0 
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THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


THEN 


(1.LE.30) THEN 


0.9454+(1I-25)*0.01/5.0 


(I.LE.40) THEN 
0.955+(1I-30)*0.011/10.0 
(I.LE.60) THEN 
0.966+(1-40)*0.012/20.0 
(12 EE.80)° THEN 
0.978+(1I-60)*0.006/20.0 


(I.LE.100) THEN 


0.984+(1-80)*0.003/20.0 


(I,LE,120) THEN 


0.987+(I-100)*0.003/20.0 


SUBROUTINE MLESHAPE(T,N,R, DEL, B, BNEW) 


C This subroutine returns a biased estimator (BNEW) for a Weibull 


ie, 
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shape parameter using the Newton-Raphson’s Method of Successive 
Approximation. The data parameters consist of an ascending ordered 
list of failure times (T), sample size (N), number of failed samples 
(R), tolerance for convergence (DEL) and an initial estimate of the 
shape parameter (B). 


LOGICAL DONE 

INTEGER N, R, I 

REAL*8 GFUNCT, GPRIME, B, BOLD, BNEW, T(*), DEL, 
* TERM1, TERM2, TERM3, SUM1, SUM2, SUM3, SUM4, STEP 


BNEW = B 
DONE = .FALSE. 


DO WHILE (.NOT. DONE) 


DONE = .TRUE. 


TERM] = FLOAT(N-R)*(T(R)**BNEW) 

TERM2 = FLOAT(N-R)*(T(R)**BNEW) *LOG(T(R) ) 

TERM3 = FLOAT(N-R)*(T(R)**BNEW) *LOG(T(R) )*LOG(T(R) ) 
SUM1 = 0.0 

SUM2 = 0.0 

SUM3 = 0.0 

SUM4 = 0.0 


DO 50 T= 1, R 
SUM1 = SUM1 + T(1I)**BNEW 
SUM2 = SUM2 + (T(1)**BNEW)*LOG(T(I)) 
SUM3 = SUM3 + (T(1)**BNEW)*LOG(T(I))*LOG(T(I)) 
SUM4 = SUM4 + LOG(T(I)) 
50 CONTINUE 


GFUNCT = (SUM2+TERM2)/(SUM1+TERM1) - (1.0/BNEW) 


* - (1.0/FLOAT(R) )*SUM4 

GPRIME = (1.0/(SUM1+TERM1)**2)*( (SUM1+TERM1)*(SUM3+TERM3) 
* - (SUM2+TERM2)**2 ) 
* + (1.0/BNEW**2) 


PRINT 60, GFUNCT,GPRIME, BNEW 
60 FORMAT(1X, 'GFUNCT =',F8.3,’ GPRIME =’,F8.3,’ BNEW =’, F8.3) 


Control magnitude of the marching step towards convergence 
as no more than 0.1. 


IF ((GFUNCT.LT.0) .AND. (GPRIME.GT.0)) THEN 
STEP = VMAX(-.1D0, (GFUNCT/GPRIME) ) 

ELSEIF ((GFUNCT.GT.0) .AND. (GPRIME.LT.0)) THEN 
STEP = VMAX(-.1D0, (GFUNCT/GPRIME) ) 

ELSE 
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STEP = VMIN(.1D0, (GFUNCT/GPRIME) ) 
ENDIF 


BOLD BNEW 
BNEW = BNEW - STEP 


Check for convergence of the MLE for the shape parameter B. 


IF (ABS(BOLD-BNEW) .GT.DEL) THEN 
DONE = .FALSE. 
ENDIF 


Avoid overflow error due to large MLE value caused by small 
GPRIME (slope) as GFUNCT approaches to near zero. 
STOP WHEN MAGNITUDE OF BNEW EXCEEDS 7. 


IF (BNEW.GT./7.0) THEN 


BNEW = BOLD 
DONE = .TRUE. 
ENDIF 
ENDDO 
RETURN 
END 
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FUNCTION CHISQD(P,N) 


Modified version of Algorithm 451 from Comunications of the ACM 
Aug 1977 Vol.16 No.8 


This function evaluates the quantile at the probability level P 
(left tail area) for the Chi-square distribution with 
N degrees of freedom. 


REAL*8 P 

REAL X 

INTEGER IF 

DIMENSION C(21), A(19) 

DATA C/ 1.565326E-3, 
1.060438E-3, 
-6.950356E-3, 
3232935 
2.277679E-2, 
-8.986007E-3, 
-1.513904E-2, 
2.530010E-3, 
Ie45Oligeage 
5 .169654E-3, 


+e ee He HF Fe 
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212153761 aoe 
1.128186E-2, 
2.607083E-2, 
-0.2237368, 
9.780499E-5, 
-8.426812E-4, 
3.125580E-3, 
-8.553069E-3, 
1.348028E-4, 
0.4713941, 
1.0000886 / 
DATA A/ 1.264616E-2, 
-1.425296E-2, 
1.400483E-2, 
-5.886090E-3, 
-1.091214E-2, 
-2.304527E-2, 
3. 1354013), 
-2.728484E-4, 
-9.699681E-3, 
13ee7 2-2. 
2.618914E-2, 
20h 22222228 
5.406674E-5, 
3.483789E-5, 
=7 274761b-4, 
3.292181E-3, 
=8 7207 IGE 
0.4714045, 
d- 
IF (N-2) 10, 20, 30 
10 CALL XFROMP(.5*(1.-P),X,IF) 
CHISQD = X 
CHISQD = CHISQD*CHISQD 
RETURN 
20 CHISQD = -2.*LOG(1.-P) 
RETURN 
30 F=N 
Fl = 1./F 
CALL XFROMP(P,X, IF) 
T =X 
F2 = SQRT(F1)*T 
IF (N.GE. (2+INT(4.*ABS(T)))) GO TO 40 
CHISQD = (((((((C(1)*F2+C(2) )*F2+C(3))*F2+C(4) ) *F2 
+C(5))*F2+C(6) )*F2+C(7))*F1+((((((C(8)+C(9) *F2) *F2 
+C(10) )*F2+C(11) )*F2+C(12) )*F2+C(13) )*F2+C(14) ))*F1+ 
(((((C(15)*F2+C(16)) *F2+C(17) )*F2+C(18) )*F2 
+C(19) )*F2+C(20) )*F2+C(21) 
GO TO 50 
40 CHISQD=(((A(1)+A(2)*F2)*F1+(((A(3)+A(4) *F2) *F2 


eee eee HH 
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+A(5))*F2+A(6)))*F14+(((((A(7)+A(8) *F2) *F2+A (9) ) *F2 
+A(10) )*F2+A(11) )*F2+A(12)))*F1+(((((A(13) *F2 
+A(14))*F2+A(15) )*F2+A(16))*F2+A(17) )*F2*F2 
+A(18))*F2+A(19) 

50 CHISQD = CHISQD*CHISQD*CHISQD*F 

RETURN 

END 


es e® ew®eew eweeweeeeeBeweeereeneteweeewe=ee®ewse®Feeeeq@eoeede@eeeeqe2e2@2cqnree2e88¢qe¢e2e2@e2eoe22e2ooB2@e8eeqe2eoeoeree2e2@2e@¢@oe2qsnse@e@#eF8#2"2FB.222.£:;»«@ ws eweeeeee = = 


% % % % 


SUBROUTINE XFROMP(P,X, IFAULT) 
Aleorithm AS 24 J.R.STAT:SOC. C. (1969) Vol.18. No.3. 


This subroutine computes the standard normal deviate X for 
the specified left tail area P. 


REAL*8 P 

DIMENSION A(5) 

DIMENSION CONNOR (17), HSTNGS(6) 
DATA CONNOR/ 8.0327350124E-1/7, 
-4483264644E-15, 
.46682/70103E-14, 
.9554295164E-13, 
.9477940136E-12, 
.3507027951E-11, 
POSIZZ2Z1037E-9, 
.3122532964E-8, 
-4503852223E-7, 
.4589169001E-6, 
Hoe 27 152208 - >: 
.0683760684E-4, 
~9/5757/5758E-4, 
~6296296296E-3, 
-3809SZ5810E-2, 
ml, 

53393553353 007 


Fe ee FF + ee eH HH HF FH HF SF 
OONnNFN RP RRP rrr wOUWN PF 


DATA RTHFPI / 1.2533141373 / 
DATA RRT2PI / 0.3989422804 / 
DATA TERMIN / 1.0E-11 / 


DATA HSTNGS / 2.515517, 
0.802853, 
010328, 
432788, 
189269, 
001308 / 


+ + + + 
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Get 


IFAULT = 1 
IF ((P.LE- OMO9MOR. (CPPGE. 1.0) G0" teas 
IFAULT = 0 


first approximation XO to deviate by Hastings’ formula 


B = P 
IF(B.GT.0.5) B= 1.0 - B 


F = - LOG(B) 
E = SQRT(F+tF) 
XO = -E + ((HSTNGS(3)*E+HSTNGS (2) )*E+HSTNGS (1) )/ 


* (( CHSTNGS (6) *E+HSTNGS (5) ) *E+HSTNGS (4) ) *E+1 .0) 


IF (XO. LT ORONMeS TOMl 
XO = 0.0 

PO = 0.5 

Xl = -RTHFPI 

GO 107 


Find the area PO corresponding to XO 


1 


(1) 


(2) 


Y = XO*%*2 
IF °(X0'. LES- 1 9)6Go TO 3 
Y = -0,5*Y 


series approximation 


PO = CONNOR(1) 
DO 2 L=2,17 
PO = PO*Y + CONNOR(L) 


PO = (PO*Y+1.0)*XO 

Xl = -(PO+RTHFPI)*EXP(-Y) 
PO = PO*XRRT2PI + 0.5 

GO TO 7 


continued fraction approximation 


Z=1.0/Y 

A(2) = 1.0 
AG) 1 6 
A(4) =OZ + 1.0 
A(5) = 1.0 
Wee 2.0 

DO 6 L=loe2 
DO 5 J=1,2 
Rows o 

ne ee 
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5 ACK) = ACKA) + A(K)*W*Z 


6 W= W+ 1.0 
APPRXU = A(2)/A(3) 
APPRXL = A(5)/A(4) 
C = APPRXU - APPRXL 
IF (C.GE.TERMIN) GO TO 4 
Xl = APPRXL/XO 
PO = -X1*RRT2PI*EXP(-0.5*Y) 


Get accurate value of deviate by Taylor Series 
(X1, X2, X3 are derivatives for the Taylor Series 


ja = F + LOG(FO) 
X2 = XO*X1*X1 -X1 
X3 = X1*¥*3 + 2.0*XO*X1*X2 -X2 
X = ((X3*D/3.0+X2)*D/2.0+X1)*D + XO 
IF (P.LE.0.5) GO TO 100 
a=" .X 
100 RETURN 
END 


FUNCTION VMAX(X,Y) 

REAL*8 X, Y 

Tho Ct. Ch. mernEN 
VMAX = X 

ELSE 
VMAX = Y 

ENDIF 

RETURN 

END 

FUNCTION IMAX(X, Y) 

INTEGER X, Y 

IF (X.GT.Y) THEN 
IMAX = X 

ELSE 
IMAX = Y 

ENDIF 

RETURN 

END 


FUNCTION VMIN(X, Y) 

REAL*8 X, Y 

IF (X.LT.Y) THEN 
VMIN = X 

ELSE 
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VMIN = 
ENDIF 
RETURN 
END 


FUNCTION 


Y 


IMIN(I,J) 


INTEGER I, J 


TE Giger 
IMIN = 
ELSE 
IMIN 
ENDIF 
RETURN 
END 


.J) THEN 
I 


SUBROUTINE GETP(N,ALPHA,NEWP) 
INTEGER N 

REAL*8 ALPHA, OLDP, NEWP, TOL 
LOGICAL DONE 

OLDP = ALPHA**(1.0/FLOAT(N) ) 


NEWP = 1 


0 


TOL = 0.0001 
DONE = .FALSE. 


DO WHILE 
GFUNCT 
GPRIME 
NEWP = 


(NOT. DONE) 

= N*OLDP**(N-1) - (N-1)*OLDP**N - ALPHA 

= N*(N-1)*OLDP**(N-2) - N*(N-1)*OLDP**(N-1) 
OLDP - (GFUNCT/GPRIME) 


IF ((ABS(NEWP-OLDP).LE.TOL) .OR. (ABS(GFUNCT) .LE.TOL)) THEN 


DONE 
ENDIF 
OLDP = 

END DO 
RETURN 
END 


= .TRUE. 


NEWP 


4. Program Output. (OUT2.DAT) 


The result 


for the simulation run based on the input parameters specified in 


IN2.DAT are computed and written to the file OUT2.DAT. A sample of this file is 


as follows. 


OUT2.DAT : Output File of the RETP2 simulation 
after 1000 replications 


COMMENTS : 8 COMPONENT IN SERIES 


DF 


= NINT (1.3 * 2 * (1 + NFC)) 
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Input Parameters: 


I TY(I) PARM1(T) 


CONDO FWHM F 


ISEED 


16807 .0 


NONNY NY Free 
ooooqoo 0 © 


NCOMP 


PARM2 (I) 
.00500 1.00000 
.00500 1.00000 
.00500 1.00000 
.00500 1.00000 
.01000 2.00000 
.01000 2.00000 
.01000 2.00000 
.01000 2.00000 


ALPHA 


8 0.2000 0.01000 


Output Parameters for the LAST Replication: 


I 


COM DOE WD 


Cut - 


: 


CONMDU LF WN FEF 


NF(T) 


Set Data: 


PRP RP RPP RP 
ONNDUPWNPE 


RS 


0.8269590 


ET(I) 


. S400000E+04 
. 5400000E+04 
. S4400000E+04 
. S400000E+04 
0237 O76 Eto 
.4825114E+06 
.25939978E+05 
.3708000E+07 


oO OO 020 0 © 


Component List 


Ooo OO © © OC @ 


So oC ca. © 2 oa 
Oo Coo 0 0 2 @ © 
oO oo oo 0 a 2 


ELMAX (M) 


0.0066667 


ooo ooo 0C Oo 


TOL NCS TCN 
8 3 
UICr) creme NC(L)  NF(I) 
5.00 5400.00 20 36 
5.00 5400.00 20 23 
5.00 5400.00 20 28 
5.00 5400.00 20 25 
15.00 2700.00 20 20 
15.00 2700.00 20 20 
15.00 2700.00 20 20 
15.00 2700.00 20 20 
ELM (I) ER(I) EBETA(1) 
0.0066667 1.0000000 0.0000000 
0.0042593 0.6388889 0.0000000 
0.0051852 0.7777778 0.0000000 
0.0046296 0.6944444 0.0000000 
0.0008407 0.1261042 1.5441284 
0.0000414 0.0062175 2.1482593 
0.0007874 0.1181112 1.5843514 
0.0000054 0.0008091 2.5186896 
RELI REL2 (M) 
0 0 O 0.975309908 0.955163479 
0 0 O 0.975309908 0.971117675 
0 0 O 0.975309908 0.964950144 
0 0 O 0.975309908 0.968645990 
0 O O 0.977751195 0.927053154 
0 O O 0.977751195 0.981007099 
0 O O 0.977751195 0.923940539 
0 0 O 0.977751195 0.993218899 
LMU (M) RSLOW LEVEL 
0.0091745 0.7754698  0.9940000 


87 


APPENDIX D : Evaluation of Subroutines and Functions 


RANDOM NUMBER GENERATOR 
(LRNDPC) 
Evaluation 


One thousand uniform random real numbers between 0 and 1 are generated 
using the random number generating routine LRNDPC. From these uniformly 
distributed numbers, 1000 exponential (with scale parameter 1) numbers and 1000 
Weibull (with scale parameter 1 and shape parameter 2) numbers were generated. 


Uniform Random Variate 


Plot of U(N+1) vs U(N) 
for 1000 Uniform RVs from LRNDPC 


U(N+1) 





0 Ot O2 OS 04 O8 O86 OF O8 OS 1 
U(N) 


Figure 1: Uniform RVs generated by LRNDPC 


Figure 1 above shows a plot of 1000 uniform real numbers against their 
predecessors. The uniformity of the distribution of points over the state space 
confirms LRNDPC’s adequacy in generating uniform random numbers. 
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Exponential and Weibull real numbers were generated using these 1000 
Uniform(0,1) random numbers. The cumulative histograms of these resultant 
random variates were compared with their respective theoretical cumulative 
distribution functions (cdfs). 


Exponential Random Variate 


Distribution Fitting 
for 1000 Exp(1) RVs by LRNDPC 


Cumulative Proportion 





0 1 2 3 4 § 6 7 8 9 10 
Value of Random Variable 


Figure 2 : Exponential RVs generated by LRNDPC 


Figure 2 shows the close distribution fit between the theoretical cdf (line) and 
the cumulative distribution of exponential RV generated using LRNDPC. 


F(t) 
mee é 


exp(-Af) 
-< IntF(e) 
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Weibull Random Variate 


Distribution Fitting 
for 1000 Wel(1,2) RVs by LRNDPC 
(le. Scale = 1, Shape = 2) 


Cumulative Proportion 





0 05 1 15 2 25 3 3.5 4 
Value of Random Variable 


Figure 3 : Weibull RVs generated by LRNDPC 


Figure 3 shows the close distribution fit between the theoretical cdf (line) and 
the cumulative distribution of Weibull RV generated using LRNDPC. 


exp{-(Ar)*} 
~~ [intFo)]* 


F(t) 


f 
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Plot of Unbiasing Factor B(N) vs N 
for Weibull Shape Parameter Estimation 


N = Test Sample Size 
B(N) = Unbiasing Factor for MLE 


Plot of BN vs N 
(Welbull Shape Parameter Estimation) 


0.95 
0.9 
0.85 


08 


Unbiasing Factor BN 


0.7 


0.65 


0 2 «0 80 100 120 


60 
Sample Size N 
Figure 1: B(N) vs N 


The function BN(N) returns the linear-interpolated values of the unbiasing 
factor for the raw MLE 8 for both RETP1 and RETP2. 


Evaluation of Subroutine CHISQD and XFROMP 
The x’ statistics for 1 to 499 degrees of freedom for a@ values of 0.1 and 0.2 are 


generated using the routines CHISQD and XFROMP. These outputs matched those 
tabulated in the mathematical tables of any general textbook on statistics. 
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I ee ce ee ee ee 


oan nwa 


10 
abel 
12 
13 
14 


15 
16 
17 
18 
19 


20 
21 
22 
23 
24 


25 
26 
27 
28 
29 


30 
31 
32 
33 
34 


35 
36 
37 
38 
39 


40 
41 
42 
43 
44 


45 
46 
47 
48 
49 


15.99 
28.41 
40.26 
5181 


63.17 
74.40 
S2-50 
96.58 
107. 57 


118.50 
129.39 
140.23 
151,05 
161.83 


172.58 
183.31 
194.02 
204.70 
21507 


226.02 
236.65 
247.27 
257.88 
268.47 


279.05 
289.62 
300.18 
310.72 
321.26 


331.79 
342.31 
352.82 
363.32 
373.82 


384.31 
394.79 
405.26 
415.73 
426.19 


436.65 
447.10 
457.54 
467.98 
478.42 


488.85 
499.27 
509.69 
520.11 
530.52 


2.71 
eee 
29.62 
41.42 
22, 99 


64.30 
Took 
86.64 
97.68 
108.66 


119-58 
130.47 
141.32 
SZ 2 
162.90 


173.66 
184.38 
195.09 
205.77 
216.44 


227.09 
237.72 
248.33 
258.94 
269.53 


280.11 
290.67 
301.23 
311278 
322.31 


332.84 
343.36 
353.87 
364.37 
374.87 


385.35 
395.84 
406.31 
416.78 
427.24 


437.69 
448.14 
458.59 
469.03 
479.46 


489.89 
500.32 
510.74 
S2u.1> 
531.56 


174.73 
1o5.45 
196.16 
206.84 
217.20 


228.15 
238.78 
249.40 
260.00 
270.59 


281.16 
291.73 
302.29 
312.83 
323.37 


333.89 
344.41 
354.92 
365.42 
375.92 


386.40 
396.88 
407.36 
417.82 
428.28 


438.74 
449.19 
459.63 
470.07 
480.51 


490.93 
501.36 
511.78 
522.19 
532.60 


176.88 
187.60 
198.29 
208.97 
219.63 


230.28 
240.90 
2512 52 
262.12 
272.70 


283.28 
293.84 
304.40 
314.94 
325.47 


336.00 
346.51 
357.02 
367.52 
378.01 


388.50 
398.98 
409.45 
419.92 
430.38 


440.83 
451.28 
461.72 
472.16 
482.59 


493.02 
503.44 
513.86 
524.28 
534.69 


252.58 
263.18 
273.76 


284 .34 
294.90 
305.45 
315.99 
326.53 


337.05 
347.56 
358.07 
368.57 
379.06 


389.55 
400.03 
410.50 
420.96 
431.42 


441.88 
452.32 
462.77 
473.20 
483.63 


494.06 
504.49 
514.90 
525.32 
S35 0rS 


125.04 
135.90 
146.72 
157.22 
168.28 


179.02 
189.74 
200.43 
Pall el 
221.76 


232.40 
243.03 
253.64 
264.24 
274.82 


285.39 
295.95 
306.51 
317.05 
327.58 


338.10 
348.62 
359.12 
369.62 
380.11 


390.60 
401.07 
411.54 
422.01 
432.47 


442.92 
453.37 
463.81 
474.25 
484.68 


495.10 
5305.26 
515.95 
526 .36 
536.77 


v2 


233.47 
244.09 
254.70 
265.29 
275.88 


286.45 
297.01 
307.56 
318.10 
328.63 


339.15 
349.67 
360.17 
370.67 
381.16 


391.64 
402.12 
412.59 
423.05 
433.51 


443.96 
454.41 
464.85 
475.29 
485.72 


496.15 
506.57 
516.99 
527.40 
5a7 260 


i272) 
138.07 
148.89 
159.67 
170.43 


16117 
191.88 
202197 
213.24 
223.89 


234.53 
245.15 
Z55.76 
266.35 
276.94 


287.51 
298.07 
308.61 
319.15 
329.68 


340.20 
350.72 
361.22 
371.72 
382.21 


392.69 
403.17 
413.64 
424.10 
434.56 


445.01 
455.46 
465.90 
476.33 
486.76 


497.19 
507.61 
518.03 
528.44 
556.05 


341.26 
351.77 
362.27 
372.77 
383.26 


393.74 
404.21 
414.68 
425.15 
435.60 


446.05 
456.50 
466.94 
477.38 
487.81 


498.23 
508.65 
519.07 
529.48 
539.89 


a5 
16 
17 
18 
19 


20 
Zi 
ZZ 
23 
24 


Z5 
26 
27 
28 
29 


30 
31 
32 
33 
34 


35 
36 
37 
38 
39 


40 
41 
42 
43 
44 


45 
46 
47 
48 
49 


3. 
22), 
36. 
WAL, 


47 


58. 
Be: 
79. 
90. 
05 


101 


ilabt 


122. 
132. 
34 
153. 


143 


164. 
174. 
185. 
195. 
206. 


Z16. 
227. 
ZS? ; 
.83 
ZG. 


247 


268. 
278. 
289. 
299. 
310. 


320. 
330. 
.07 
.40 
361. 


341 
351 


372. 
382. 
392. 
.99 
413. 


402 


423. 
433. 
444, 
47 
464, 


454 


475. 
485. 
495. 
505. 
1G. 


44 
04 
25 


16 
97 
71 
41 


se7 


25 
81 


85 


35 
83 
29 
74 
18 


61 
03 
43 


22 


60 
97 
34 
70 
05 


40 
74 


73 


05 
37 
68 


29 


59 
89 
18 


7s) 


03 
31 
59 
86 
13 


59. 
70), 
80. 
1. 
OZ: 


2s 
12d 
133. 
144, 
154. 


165. 
175. 
186. 
196. 
<2o 


207 


217 


269. 
280. 
290. 
300. 
oui 


SZ. 
331. 
342. 
SoZ. 
362. 


373. 
383. 
393. 
.02 


404 


414. 


424, 
434. 
445. 
455. 
465. 


476. 
486. 
496. 
506. 
S17. 


25 
05 
79 
47 
12 


73 
31 
86 
39 
90 


40 
88 
34 
79 


.65 
228. 
238. 
248. 
259. 


07 
47 
87 
26 


64 
01 
37 
73 
09 


43 
Ta 
a 
44 
76 


08 
40 
rat 


32 


62 
91 
21 
50 
78 


06 
34 
62 
89 
16 


60. 
Lae 
81. 
92. 
103. 


113. 
124. 
134. 
145. 
155: 


166. 
176. 
167. 
83 
208. 


197 


218; 
229. 
239. 
249. 
.29 


260 


270. 
Zei- 
291. 
301. 
312. 


322. 
332. 
343. 
353. 
363. 


374 


425. 
435. 
446. 
456. 
466. 


477 
487 
497 
507 
518 


33 
13 
86 
54 
18 


79 
36 
91 
44 
95 


45 
92 
38 


27 


69 
11 
51 
91 


67 
05 
41 
a7 
12 


47 
81 
14 
47 
79 


oat 
384. 
394. 
405. 
415. 


43 
74 
05 
35 


65 
94 
24 
52 
81 


.09 
no? 
64 
.92 
.19 


GL 
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114. 
Zor 
135. 
146. 
7. 


167, 
177. 
188. 
198. 
acca 


209 


219. 
230. 
240. 
250. 
261. 


271 
282 
292 


323. 
333. 
344, 
. 50 
.83 


354 
364 


375. 
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APPENDIX E: Tabulated Run Results for RETP1 
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Table 1A : 8 Exp in Series, RS = 0.931 (Hi) 
min = 0.0002 f/hr, max = 0.0016 f/hr, UT = 10 hrs 


Test Deg of Measures of Accuracy 
S/N Plan Freedom 
RSLOW LEVEL 


2°NFC 
(80) 
0.2 0.919 0.9 


0.1 0.906 1.000 


0.955 


Pos | 05s | oe 
Tos [osm | ane 


0.932 0.747 


2*NFC 0.927 0.955 
Coz | ons | oni 























Test 5 until 
5 failed. 


NFC=40 


Test 15 until 
15 failed. 


NFC= 120 


Test 15 until 
11 failed. 


NFC=88 0.921 0.916 


a 
a 


w)5) 


Table 1A : 8 Exp in Series, RS = 0.931 (Hi) (Cont...) 
min X = 0.0002 f/hr, max } = 0.0016 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
RSLOW _ LEVEL 


Test 15 until 2*NFC 
fe a% 
NrC=s 


Test 15 until 
3 failed. 
Poa 0912 0.975 


NFC=24 
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Table 1B: 8 Exp in Series, RS = 0.803 (Lo) 
min X = 0.001 f/hr, max } = 0.0045 f/hr, UT = 10 hrs 


Deg of 
Freedom 
{|  RSLOW | | _LEVEL 


———— ee 


2NFC 077 
(80) 
a2 | 079s 
pos | ome | 1.000 


0.1 : 
a 
or | 07 | ose 
Pee | eon [are 












Test 15 until 
15 failed. 









NFC= 120 





11 failed. 


NFC=88 





oF 


Table 1B: 8 Exp in Series, RS = 0.803 (Lo) (Cont...) 
min A = 0.001 f/hr, max X = 0.0045 f/hr, UT = 10 hrs 


Deg of 
Freedom 
RSLOW LEVEL 


Test 15 until 
7 failed. 


NFC=56 


Test 15 until 
3 failed. 


NFC=24 


2*(NFC- 
NCOMP) 
32 
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Table 2A: 8 Wei in Series, RS = 0.980 (Hi) 
min = 0.001 f/hr, max ) = 0.008 f/hr, UT = 10 hrs 


Test Deg of Measures of Accuracy 
S/N Plan Freedom 
RSLOW 
2°NFC os | ose 
(80) 
0.2 0.989 









Test 5 until 
5 failed. 








NFC=40 0.937 0.994 


roa econ 
[oz [os [ese 
2*NFC 
a 


ee 


11 failed. (176) 
0.2 0.977 0.860 








Test 15 untl 
15 failed. 









NFC= 120 





NFC=88 0.973 0.914 


re [ oo [oe 
ore | ox | oom | cam 


NCOMP) 
160) 0.2 0.979 0.819 





o 


Table 2A: 8 Wei in Series, RS = 0.980 (Hi) (Cont...) 
min = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
_— | RSLOW | LEVEL LEVEL 


Test 15 until 0.987 
7 failed. 
0981 0.779 


NFC=56 0.985 0.839 


Test 15 until 
3 failed. 


NFC=24 
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Table 2B : 8 Wei in Series, RS = 0.832 (Lo) 
min X = 0.003 f/hr, max A = 0.024 f/hr, UT = 10 hrs 


oe 2s 
a 
Cor [ome | mn 
rez | ooo 
0.876 


0.825 




























Test 5 until 
5 failed. 


NFC=40 


Test 15 until 
15 failed. 


NFC=120 


NS 

5, 

-) 
S 


Test 15 until 
11 failed. 


NFC=88 


101 


Table 2B : 8 Wei in Series, RS = 0.832 (Lo) (Cont...) 
min X = 0.003 f/hr, max X = 0.024 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
_RSLOW LEVEL 


Test 15 until 2*NFC 
0.2 0.857 
Nos 


0.2 0.888 


nrc rex [eas | or 


ee 
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Table 3A: 4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) 
min X = 0.002 f/hr, max X = 0.008 f/hr, UT = 10 hrs 


Test Deg of 
S/N Plan Freedom 
RSLOW _ LEVEL 
2*NFC 
(80) 
ie 0.978 0.905 
| ssh 
eae ee 


ee 
863 


0 










Test 5 until 
5 failed. 


= 


NFC=40 


Test 15 until 
15 failed. 


NFC=120 0.98% 0.971 


0.1 
Toe | ene | ome 


el ss 


She Ml Test 15 until 
11 failed. 


NFC=88 0.979 0.951 


ia ea 


0.982 0.591 
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Table 3A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) (Cont...) 
min X = 0.002 f/hr, max 4 = 0.008 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 


Test 15 until 
3 failed. 


NFC=24 
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Table 3B : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.809 (Lo) 


min A = 0.002 f/hr, max X = 0.008 f/hr, UT = 10 hrs 
1 
















Deg of Measures of Accuracy 
Freedom 


RSLOW 


2*NFC 0.788 
0.754 
0.794 
0.805 
0.823 0.838 


1 
0.2 0.817 0.787 


3 | Test 15 until 2*NFC 0811 
11 failed. (176) 3 
0.807 


REVEL 





Test 5 until 
5 failed. 


0.961 


NFC=40 


* | 


0.987 
0.910 
0.874 


Test 15 until 
15 failed. 


NFC= 120 


NFC=88 0.1 0.797 0.962 


oars 
menee 


NCOMP) 
(160) 0.2 0.822 0.647 
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Table 3B : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.809 (Lo) (Cont...) 
min A = 0.002 f/hr, max = 0.008 f/hr, UT = 10 hrs 


Deg of Measures of Accuracy 
Freedom 
_RSLOW {| LEVEL ~ 
Test 15 until 2*NFC 0.814 0.872 
7 failed. | (112) 
oa "er 
NFC=56 4 0.792 0.963 
hee 5081 


0.550 


Test 15 until 
3 failed. | 
0.755 _ 


0.780 
fe ah 
a 
ae 
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APPENDIX F : Tabulated Run Results for RETP2 
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Table 4A : 8 Exp in Series, RS 
= 0.001 f/hr, UT 


0.961 (Hi) 
5 hrs 


Degrees. | K / E[NFC] 
of (TT) 
Freedom RSLOW _ LEVEL 


2(L+NFC) | 025/12 0851 
225 
aaa 08st 


md Fe a 
a pee ee 
on 


30.0 / 144 0.906 
(27000) han r 
fe UL2 0.961 0.804 
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Table 4A : 8 Exp in Series, RS = 0.961 (Hi) (Cont...) 
d = 0.001 f/hr, UT = 5 hrs 


K / E[NFC] Measures of Accuracy 
aad | rstow | LEVEL 
Freedom _ RSLOW —_ LEVEL 


13 0.25 / 12 0.950 
— me [on | os oa 
0.2 0.935 0.851 
oe [Sar 
450 
_ 0857 


10 / 48 0.51 
(900) 
0.935 0.941 


2.0 / 9.6 
(1800) 


0.989 
0.948 


3.0 / 14.4 
(2700) 


0 
0 
02 | 090 
0 0.949 0,969 
oy oe Z 
(3600) 

ee 
5.0 / 24 
(4500) 

0.985 
10.0 / 48 
(9000) 

ls 
20.0 / 96 0.949 1. 


000 
(18000) 
0.2 0.950 1.000 
30.0 / 144 0.950 


(27000) 
0.2 0.949 1.000 | 


2 

1 
2 
1 
2 
0.1 
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Table 4B : 8 Exp in Series, RS = 0.819 (Lo) 
= 0.005 f/hr, UT = 5 hrs 


K / EINFC 
we | vow] even 
eS RSLOW | LEVEL 
2*(1+NFC) | 025/12 0.851 
| 45 ——- 
—— asst 
05 / 24 0887 
(90) ope 


si 
20 / 96 
ae 
5.0 / 24 
ap 
0 


ae 0.818 
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Table 4B : Ome XE I SCUIES, ISO 
X = 0.005 f/hr, UT 


K / EINFC) 
(TT) 
| RSLOw | _LEVEL 


| (eee 
ee x 
5.0 / 24 
<i 
ene Ean 0s asp 
Cc) or2 0.772 1.000 


0.819 (Lo) (Cont...) 
5 hrs 


























J 


Table 5A: 8 Wei in Series, RS = 0.956 (Hi) (*) 
= 0.005 f/hr, UT = 15 hrs 


Degree eee 
eee ai 
2*(+NFO) 
ull 

Toa | oom | cae 


1.0 / 48 
(180) 


2.0 / 9.6 


360 
: 0.952 0.854 





(540) 


4.0 / 19.2 
(720) 


5.0 / 24 
(900) 


10.0 / 48 
(1800) 


20.0 / 96 
(3600) 


30.0 / 144 
(5400) 





(*) 20 test items for each Weibull component. 
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Table SA : 8 Wei in Series, RS = 0.956 (Hi) (*) (Cont...) 
= 0.005 f/hr, UT = 15 hrs 


of (TT) 
Freedom |  ———ssdTS——idT_~C«éRSLOW RSLOW _ LEVEL — 
13° 025 / 12 0258 
2*(1+NFC) (45) 
0.2 0.999 0.224 
(i595 02.4 0.983 0.635 
(90) 
0.2 0.983 0.593 
190 / 48 0956 
(180) 
0.2 0.930 0.866 
20/96 
(360) 
0.2 0.939 | 0.956 
(540) 


(720) 

900 

oe | os | ae 
100 / 48 0.981 
1800 

= 0.981 
20.0 / 96 0.981 
3600 

——— 0.981 


(5400) 
2 0.930 0.981 








3 


Table SB: 8 Wei in Series, RS = 0.835 (Lo) (*) 
= 0.01 f/hr, UT = 15 hrs 


Degrees SGT |g; Meeaettee 
= _ 
marry | 0/2 | on | ome | oz | 
= 

05 / 24 

= 

10 / 48 

‘r 

20 / 96 

2% 

30 / 144 

sojwa | ox | om | ome 

50 / 2 

i 

100 / 48 

_ 

20.0 / 96 

a 

30.0 / 144 


(2700) 
0.813 0.908 





114 


Table SB: 8 Wei in Series, RS = 0.835 (Lo) (*) 
X = 0.01 f/hr, UT = 15 hrs 


K / EINFC 
TT 
Freedom = RSLOW LEVEL 
025 / 12 0957 
2*(1+NFC) (22.5) 
02 | 095s 
05 / 24 0572 
45 
z. 
1.0 / 48 
90 
¥ 
20 / 96 


180 
(270) 
(360) 
(450) 

0.2 0.766 0.986 
(900) 

0.2 0.766 0.986 
20.0 / 96 0.784 0.989 
(1800) 

0.2 0.766 0.986 


30.0 / 144 0.989 
(2700) | 





0.766 0.986 


NS) 


Table 6A: 


4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) 


(exp) = 0.001 f/hr, UT(exp) = 5 hrs 


(wei) = 


Degrees K / E{[NFC] 
of 


Freedom 


2*(1+NFC) 





3.0 / 14.4 
(2700) 
(540) 


4.0 / 19.2 
(3600) 
(720) 


20.0 / 96 
(18000) 
(3600) 


30.0 / 144 
(27000) 
5400 
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0.005 f/hr, UT(wei) = 


15 hrs 


RSLOW LEVEL 


0.956 0.862 





Table 6A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) (Cont...) 
(exp) = 0.001 f/hr, UT(exp) = 5 hrs 
(wei) = 0.005 f/hr, UT(wei) = 15 hrs 


me K / E[NFC] Measures of Accuracy 
Tris eve 
w|reedom | tive) | | stow | Lever 
135 
2*(1+NFC) 


3.0 / 14.4 
(2700) 
(540) 
4.0 / 19.2 : 0.948 
(3600) 

0.943 0.965 
(18000) 





ely 


Table 6B : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.827 (Lo) (*) 
(exp) = 0.005 f/hr, UT(exp) = 5 hrs 
(wei) = 0.010 f/hr, UT(wei) = 15 hrs 


i K / E{[NFC] Measures of Accuracy 
TT (exp) 
RSLOW LEVEL 


Bs" TT (wet) 


2*(1+NFC) 


a) on 
ce 
(540) 
(720) 
0.822 


: 0.902 
igs [ae 


(5400) 
2700 0.821 | Se ORS4 ee 854 
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Table 6B : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.827 (Lo) (*) (Cont...) 
(exp) = 0.005 f/hr, UT(exp) = 5 hrs 
(wei) = 0.010 f/hr, UT(wei) = 15 hrs 


ca K / E[NFC] Measures of Accuracy 
TTve | RsLow | LEVEL 
| ee 


TS 
2*(1+NFC) 


= — 


(540) 
4.0 / 19.2 0.792 0.939 
(720) 


Fe 


(5400) 
(2700 Sosa 775 





iy, 
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